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Abstract 

The simplest high gain free electron laser (FEL) amplifier concept is proposed. 
A klystron amplifier has the useful property that the various electronic processes 
take place in separate portions of the amplifier, rather than overlapping as in FEL 
amplifier with an uniform undulator. The klystron consists of two fundamental 
parts: succession of 2-3 cascades (modulator), and an output undulator (radiator) 
in which the modulated electron beam coherently radiates. Each cascade consists of 
uniform undulator and dispersion section. Unlike distributed optical klystrons, we 
have a high gain per cascade pass. This has a few consequences. First, klystron gain 
does not depend on the bunch compression in the injector linac, i.e. maximum gain 
per cascade pass at high peak beam current is the same as at low peak beam current, 
without compression. Conventional, short- wavelength FEL amplifier and distributed 
optical klystron require electron beam peak current of a few kA. Second, there are 
no requirements on the alignment of the cascade undulators and dispersion sections, 
because, in our (high gain) case, there is no need for radiation phase matching. There 
are applications, like XFELs, where the unique properties of high gain klystron FEL 
amplifier are very desirable. Such a scheme allows one to decrease the total length of 
the magnetic system. In this paper we discuss also implementation of the proposed 
SASE FEL scheme for the frequency multiplier, for femtosecond experiments, and 
multi-user soft X-ray facility. 
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1 Introduction 



High gain FEL amplifiers are of interest for a variety of potential applications 
that range from X-ray lasers [1,2] to ultraviolet MW-scale industrial lasers 
[3]. There are various versions of the high gain FEL amplifier. A number of 
high gain FEL amplifier concepts may prove useful for XFEL applications. 
Two especially noteworthy ones are the FEL amplifier with a single uniform 
undulator [1,2] (see Fig. 1) and the distributed optical klystron [4-7] (see 
Fig. 2). The high gain cascade klystron amplifier described in this paper is 
an attractive alternative to other configurations for operation in the X-ray 
wavelength range (see Fig. 3). 

Let us first investigate qualitatively the processes that take place in any FEL 
amplifier; this investigation should also enable us to distinguish somewhat 
more accurately between the so-called distributed optical klystron concept and 
the new high gain klystron amplifier concept. In any FEL, several processes 
must take place, either separately or simultaneously. The first process that 
must occur in any FEL is energy modulation. Energy modulation means the 
application of radiation fields to the electrons as they move in their trajectories 
under the action of undulator magnetic fields, in such a way that the energies 
of these electrons are varied in some periodic fashion. It is then necessary to 
take advantage of this phase-dependent motion. 

The second process, then, that must occure in a FEL is the conversion of 
this variation in the motion of the electrons into a usable form. It is neces- 
sary to turn the variations in the electron energy into a radiation frequency 
current, which can deliver energy to the electromagnetic field; i.e., the varia- 
tions in velocity must be converted into density variations. The ways in which 
this conversion is brought about are numerous, and FEL amplifier types are 
distinguished by their conversion mechanism. 

It may be useful to illustrate processes at this point to see how energy modu- 
lation and conversion into density modulation take place in some of the types 
of FEL amplifier. In a conventional high gain FEL amplifier with a single uni- 
form undulator energy modulation and bunching occurs both simultaneously 
and continuously along the undulator. Energy modulation occurring at any 
one point continues to affect bunching at points beyond, with additional mod- 
ulation being added as the beam moves along. The problem of electromagnetic 
wave amplification in the uniform undulator refers to a class of self-consistent 
problem. To solve the problem, the field equations and equations of particle 
motion should be solved simultaneously. When the FEL amplifier operates in 
linear regime the transverse distribution of the radiation field remains fixed, 
while its amplitude grows exponentially with undulator length. 
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Fig. 1. Schematic diagram of conventional SASE FEL 




Fig. 2. Schematic diagram of distributed optical klystron 

Perhaps the simplest FEL amplifier, from pedagogical point of view, is the high 
gain klystron amplifier. This simplicity arises from the fact that the important 
processes of energy modulation of the beam and subsequent conversion of the 
energy modulation to density modulation occur in distinctly separate regions 
of the FEL and thus may be considered in a sequential manner. This separabil- 
ity of functions aids materially in developing an understanding of modulation 
and bunching phenomena. 



In its simplest configuration the klystron amplifier consists of an input un- 
dulator (energy modulator), and an output undulator (radiator) separated 
by a dispersion section. In the input undulator an incident electron beam is 
subjected to a radiation field, which produces an energy modulation, a modu- 
lation which after dispersion section, is converted into an density modulation 
in the beam. In the output undulator the bunched beam radiates coherently. 
This process, of conversion of energy modulation to bunching by the action of 
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Fig. 3. Schematic diagram of high gain cascade klystron amplifier 

a dispersion section, is one of the characteristics of the klystron and is one of 
the reasons why this device is so useful in many connections. 

We have considered the two-undulator klystron amplifier. In some experimen- 
tal situations this simplest configuration is not optimal. For example, this 
study has shown that in the soft X-ray wavelength range the maximal (in- 
tensity) gain of this amplifier does not exceed 30-40 dB. The obvious and 
technically possible solution of the problem of gain increase might be to use a 
klystron with three or more undulators. Suppose that we have a three undu- 
lator (or so-called two-cascade) klystron, and that we introduce a very small 
signal into the first undulator. This will result in a small value of the ampli- 
tude of energy modulation, and a resulting rather small power output from 
the second undulator. Nevertheless, the radiation field amplitude in the second 
undulator will still be much greater than that in the first undulator, on ac- 
count of the considerable amplification of the first cascade. We now disregard 
the power produced by the second undulator, but consider its effect on the 
electron beam passing through its magnetic system. The considerable radia- 
tion field in the second undulator will produce a further energy modulation of 
the beam and, when it emerges on the far side of the second undulator, it will 
have a large energy modulation, which, with a second dispersion section, will 
result in optimum bunching. This bunched beam then enters the third (radia- 
tor) undulator, and produces a large amount of power. In other words, we have 
essentially combined two stages of amplification by incorporating essentially 
two klystrons in a single unit. The final power is no greater than could be 
produced by a two-undulator klystron with the same output undulator, but 
we can secure this power with a much weaker input signal. 

Let us consider now the problem of a periodically spaced sequence of undu- 
lators and dispersion sections, with the electron beam traversing each in suc- 
cession. The larger the number of cascades the longer is the amplifier length. 
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Suppose that the increase of total length (gain) of amplifier is tempered by an 
accompanying decrease in the length (gain) of cascade. From this approach at 
very large number of cascades and small gain per cascade pass we are led to 
the other device, which was the subject of the papers [4-7]: the "distributed 
optical klystron" . At first glance this is essentially the extension klystron am- 
plifier, extended to a very large number of cascades. However, as noted earlier, 
this amplifier is considerably more complicated in principle than a klystron 
amplifier with high gain per cascade pass. When the gain per cascade pass 
is small, the cascades cannot be considered independent. It is important to 
realize that the cascades are connected optically, so that the radiation pro- 
duced, for example, into the first undulator passes on into the second and 
produces additional energy modulation of the beam. It is clear that we cannot 
conveniently treat the theory of distributed klystron by the same method we 
used in the klystron case. Because of the low gain per cascade pass used for 
the distributed klystron, it cannot be approximated as scheme with separate 
functions, and the klystron amplifier theory has all based on that assumption. 
A quite different approach to the problem, which proves to be more suited 
to mathematical development of distributed klystron theory, is more along 
the lines of an analogy with the theory of conventional FEL amplifier with 
uniform undulator. The problem of electromagnetic wave amplification in the 
distributed klystron refers to a class of self-consistent problems. For instance, 
when the distributed klystron operates in linear regime its amplitude gain 
grows exponentially with cascade numbers. From a study of the rate of ex- 
ponential build-up of the electromagnetic wave we can find the gain of the 
amplifier. 

Our theory of klystron amplifier is based on the assumption that beam density 
modulation does not appreciably change as the beam propagates through the 
undulator. This approximation means that only the contributions to the radi- 
ation field arising from the initial density modulation are taken into account, 
and not those from the induced bunching. At very long klystron undulators we 
are led to device, which is the subject of the papers [8,9]. Two-segment high 
gain FEL amplifier is variant of the standard high-gain FEL amplifier with 
uniform undulator. Because of the high gain per undulator pass used for this 
modification, it cannot be approximated as scheme with separate functions. 

The klystron amplifier analysis is not only relevant from a pedagogical point 
of view. There are reasons to expect that the new concept of FEL amplifier 
discussed here can be of interest for short wavelength applications such as X- 
ray SASE FELs. Generation of X-ray radiation using klystron amplifiers has 
many advantages, primarily because of the separation of electronic processes, 
gain, and potentially low sensitivity to alignment. In principle, the gain of the 
klystron amplifier can be increased indefinitely by adding additional cascades 
along the beam. In practice, with three-cascade klystron amplifier, a gain in 
excess of 80 dB may be obtained in the soft X-ray wavelength range. At this 
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level, the shot noise of the beam is amplified up to saturation and a cascade 
klystron amplifier operates as soft X-ray SASE FEL. 

Now our intention is to develop certain physical ideas about the operation of 
the klystron amplifier, the characteristics of its radiation, and those features 
that set it apart from other FEL sources in the short wavelength region of 
the spectrum. The resonance properties of the radiation of the electron beam 
in the klystron can be understood as follows. Since a dispersion section is 
a nonresonant device, the klystron bandwidth is restricted by the resonance 
properties of the cascade undulator, that is, by its number of periods N w . 
The typical amplification bandwidth of the klystron in the soft X-ray wave- 
length range is on the order of one percent (N w ~ 10 2 ). Since we study the 
start-up from shot noise, we assume the input current to have a homogeneous 
spectral distribution. At large number of undulator periods the spectrum of 
transversely coherent fraction of radiation is concentrated within the narrow 
band, AA/A ~ N~\ 

The actual physical picture of start-up from noise should take into account 
that the fluctuations of current density in the electron beam are uncorrelated 
not only in time but in space too. Thus a large number of transverse radia- 
tion modes are excited when the electron beam enters the undulator. These 
radiation modes have different gain. Obviously, as cascade number progresses, 
the high gain modes will predominate more and more and we can regard 
the klystron as a filter, in the sense that it filters from arbitrary radiation 
field those components corresponding to the high gain modes. If we consider 
the undulator radiation from the point of view of paraxial optics, then we 
immediately see that the high gain modes are associated with radiation prop- 
agating along the axis of the undulator, as opposed to radiation propagating 
at an angle to the axis which has a low gain. Consider an electron beam of 
radius r b moving in the (cascade) undulator with length L w . The parame- 
ter Np = rg/(AL w ) can be referred to as the Fresnel number. Here A is the 
radiation wavelength. The requirement for the mode selection to be quickly 
holds at small value of the Fresnel number. Hence, for thin electron beam with 
Fresnel number of Np ~ 1, the emission will emerge in a single (fundamental) 
transverse mode and the degree of transverse coherence of output radiation 
will approach unity. 

As mentioned above, amplification in each cascade is a two-step process. Since 
we study the start-up from shot noise there is no electromagnetic field at the 
undulator entrance and the modulation of the beam density serves as the input 
signal for the klystron amplifier. The radiation field produced an ponderomo- 
tive potential across the undulator; as is now known from general FEL theory, 
this voltage modulates the energy of the electron beam. Then the energy mod- 
ulation causes modulation of the beam density in the dispersion section. If we 
can neglect beam bunching in the cascade undulator, the amplitude of en- 
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ergy modulation at the undulator exit, SS, is proportional to the resonance 
Fourier component of beam current, 5I[ n = a [n I . Here a- in is the initial beam 
microbunching at the resonant wavelength, and 1$ is the peak beam current. 
We now examine the behaviour of the electrons as they travel in a dispersion 
section. Like the situation with usual microwave klystron, the Bessel function 
factor J\(X) represents the microbunching at the dispersion section exit. The 
parameter X in our case is proportional to the product of the amplitude of 
energy modulation 5£ and the compaction factor of the dispersion section 
R 5 q. In the high gain linear regime X is much smaller then unity and we make 
the approximation Ji(X) ~ X/2. Thus, microbunching at the cascade exit is 
proportional to input signal amplitude, a out ~ R^Ioa^. When a cascade op- 
erates in the high gain linear regime the gain per cascade pass is independent 
of input signal amplitude and proportional to the product of the compaction 
factor and the peak beam current, G = a out /a in ~ R$qIq. 

If we build a klystron, we want to have small number of cascades and high 
gain. One might think that all we have to do is to get enough gain per cascade 
pass - we can always increase compaction factor i? 56 and we can always in- 
crease G, and there is no reason why we cannot use the simple two-undulator 
klystron scheme. So the limitations of a gain per cascade pass are not that it 
is impossible to build a dispersion section that has a very large compaction 
factor. In fact, one of the main problems in klystron amplifier operation is 
preventing spread of microbunching due to energy spread into the electron 
beam. This energy spread is unimportant when the compaction factor of dis- 
persion section is relatively small, but we shell see that energy spread is of 
fundamental importance in the case of very large compaction factor. To get 
a rough idea of the spread of the microbunching, the position of the par- 
ticles within the electron bunch at the dispersion section exit has a spread 
which is equal to Az ~ R^A^/j, where A7/7 is the local energy spread into 
the electron beam. And how big is Az in the spread-out microbunching? We 
know that uncertainty in the phase of the particles is about Aip ~ 2-kAz/X. 
Therefore, the largest compaction factor that we can have is approximately 
i? 56 ~ A7/ (27rA7). This is what we wanted to deduce - that the gain maximum 
per cascade pass varies as ratio G max ~ J0/A7. It is of course desirable that 
the klystron gain be maximum. An experimenter can easily tune the gain of 
the klystron by tuning the magnetic fields in the dispersion sections. It should 
be pointed out that the local energy spread plays a very important role in the 
operation of the FEL klystron amplifier, and this characteristic distinguishes 
the klystron from conventional FEL amplifier. 

Electron bunches with very small transverse emittance and high peak current 
are needed for the operation of conventional XFELs. This is achieved using 
a two-step strategy: first generate beams with small transverse emittance us- 
ing an RF photocathode and, second, apply longitudinal compression at high 
energy using a magnetic chicane. Although simple in first-order theory, the 
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physics of bunch compression becomes very challenging if collective effects 
like space charge forces and coherent synchrotron radiation forces (CSR) are 
taken into account. Wakefields of bunches as short as 0.1 mm have never been 
measured and are challenging to predict. 

The situation is quite different for klystron amplifier scheme described in our 
paper. An distinguishing feature of the klystron amplifier is the absence of ap- 
parent limitations which would prevent operation without bunch compression 
in the injector linac. According to our discussion above, the gain maximum 
per cascade pass is proportional to the peak current and inversely proportional 
to the energy spread of the beam. Since the bunch length and energy spread 
are related to each other through Liouville's theorem, the peak current and 
energy spread cannot vary independently of each other in the injector linac. 
To extent that local energy spead is proportional to the peak current, which 
is usually the case for bunch compression, the gain will be independent of 
the actual peak current. We see, therefore, that klystron gain in linear regime 
depends only on the actual photoinjector parameters. This incipient propor- 
tionality between gain and Iq/ A*y is a temptation, in designing an XFEL, to 
go to very high values of Iq/Aj and very long values of bunch length. 

Let us present a specific numerical example for the case of a 4th generation 
soft X-ray light source based on the use RF photoinjector and superconducting 
linac with duty factor 1 %. The average brilliance of a klystron facility operat- 
ing without bunch compression in the injector linac surpasses the spontaneous 
undulator radiation from 3rd generation synchrotron radiation facilities by 4 
or more orders of magnitude. Decreasing the peak current also decreases the 
peak brilliance of the SASE FEL radiation by about of factor 100, but this 
is still 6 order of magnitude higher than that of 3rd generation synchrotron 
radiation sources. 

Our studies have shown that the soft X-ray cascade klystron holds great 
promise as a source of radiation for generating high power single femtosec- 
ond pulses. The obvious temporal limitation of the visible pump/X-ray probe 
technique is the duration of the X-ray probe. At a klystron facility operating 
without bunch compression, the X-ray pulse duration is about 10 ps. This is 
longer than the timescale of many interesting physical phenomena. The new 
principle of pump-probe techniques described in section 6 offers a way around 
this difficulty Section 6 also deals with the design strategy for the multi-user 
distribution system for an X-ray laboratory. An X-ray laboratory should serve 
several, may be up to ten experimental stations which can be operated inde- 
pendently according to the needs of the user community. On the other hand, 
the prefered layout of a conventional SASE FEL is a linear arrangement in 
which the injector , accelerator, bunch compressors and undulators are nearly 
collinear, and in which the electron beam does not change direction between 
accelerator and undulators. The situation is quite different for the klystron 
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amplifier scheme proposed in our paper. Since it operates without bunch com- 
pression in the injector linac, the problem of emittance dilution in the bending 
magnets does not exist. An electron beam distribution system based on un- 
bundled electron beam can provide efficient ways to generate a multi-user 
facility - very similar to present day synchrotron radiation facilities. 

It may be wondered, after hearing of all of the wonderful properties of the 
high gain klystron amplifier, why this advance in XFEL techniques occurred 
only this year. It should be note that the XFEL which was proposed in 1982 
[10] has had nearly 20 years of development. One of the explanations is as 
follows. For years we were led to believe that local energy spread into electron 
beam which delivered from a photoinjector is relatively large. Recent analysis 
of experimental results obtained at TTF SASE FEL [11] shows that the value 
of the local energy spread should be revised. After bunch compression, it is ex- 
pected to be about 0.2 MeV (for I ~ 2.5 kA) or smaller which is significantly 
less than the previously projected value of 1 MeV. This decrease significantly 
improves operation of klystron amplifier and extends the safety margin for 
X-ray klystron facility operation. 

This paper is intended primarily for the experimental scientist or XFEL pro- 
gram manager who has to design XFELs and who wishes to undestand the 
operating principles of XFELs. The material to be discussed provides the nec- 
essary background information required for making strategic XFEL design 
decisions. In this paper the electronic processes that occur in an XFEL will 
be discussed in elementary fashion. The XFEL theory has a, not entirely un- 
deserved, reputation for difficulty and obscurity. This is due partly to the 
complexity of systems in which (shot) noise is input signal and partly to the 
unfamiliar nature of some of the mathematical tools involved (integrodifferen- 
tial equations). Klystron amplifiers are one of the simplest examples of high 
gain FEL amplifiers and therefore we shall use them as an introduction to 
the techniques of XFEL theory. In the klystron amplifier case little demand 
is made on the reader's mathematical ability. The mathematics that is in- 
volved is particularly simple, involving very simple differential equations and 
algebraic operations and no integrodifferential equations. Such introduction to 
the operating principles of XFELs has never been done before, to our knowl- 
edge. 



2 Interaction between electrons and fields in an undulator 

Klystron amplifiers amplify the input signal. To describe the klystron amplifier 
operation, we first define the initial conditions. Important practical kind of 
initial condition refers to the case when there is no radiation field at the 
undulator entrance and the modulation of the beam current density serves as 
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the input signal for the amplifier. As stated in the introductory section, in 
undulator two physical processes take place simultaneously: 

1) The modulated electron beam excites electromagnetic waves which are 
propagated along the undulator; 

2) The electric field components of these waves produce a modulation of the 
energy of the electrons. 

The approach discussed in this section consists of deriving for these two process 
separate expressions which, when joined, give quantitative information about 
the composite phenomena taking place in the undulator. We assume that the 
density modulation effect is negligible. Hence, electrons passing through the 
first undulator receive only energy modulation, but no density modulation. 
Density modulation present when the electrons enter the second undulator, 
must have developed while the electrons were traveling in the dispersion sec- 
tion. This effect will be discussed later. 



2. 1 Collective fields produced by the modulated electron beam in an undulator 

The simplest form of bunched beam is found if we have swarm of electrons, 
all moving in the same direction (say the z-direction) with the same velocity 
(say v), but in which the number of electrons per unit volume (say n) depends 
on z. Then if the charge on the electron (— e) we can easily find the current 
density at any value of z. The quantity n must really be a function of z/v — t, 
on account of the velocity v of the electrons. Thus, the current density at point 
z, at time t, is —ev z n(z/v — t), where n{z/v — t) indicates the function n of 
the argument z/v — t. Unless n is a constant independent of its argument, this 
will lead to a current density which varies with time. The simplest variation of 
n for our present purposes is a superposition of a constant, or d-c component, 
and a periodic function of z/v — t, with an angular frequency u. First we 
consider the case where n is simply a constant plus a single cosine term: 
n(z/v — t) = rio[l + a- m cosou(z/v — t)]. We must note one important fact: n 
can never become negative, and hence a m cannot be greater than unity. 

To understand the basic principles of klystron amplifier operation, let us con- 
sider the helical undulator. The magnetic field on the axis of the helical un- 
dulator is given by 

— * 

H w = e x H w cos(k w z) — e y H w sm(k w z) , 

where k w = 2it/\ w is the undulator wavenumber and e XtV are unit vectors 
directed along the x and y axes of the Cartesian coordinate system (x,y,z). 
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The Lorentz force F = —e(v x H w )/c is used to derive the equations of motion 
of electrons with charge (— e) and mass m e in the presence of the magnetic 
field. The explicit expression for the electron velocity in the field of the helical 
undulator has the form: 

v±(z) = c9 w [e x cos(k w z) — e y sin(k w z)] , 

which means that the electron in the undulator moves along a constrained 
helical trajectory parallel to the z axis. The angle of rotation is given by the 
relation 9 W = K/j = \ w eH w /(27rm e c 2 'j), where 7 = (1 — t> 2 /c 2 ) -1 / 2 is the 
relativistic factor and v 2 = v 2 + v 2 + v 2 . As a rule, the electron rotation angle 
9 W is small and the longitudinal electron velocity v z is close to the velocity of 
light, v z ~ c. 

Let us consider a bunched relativistic electron beam moving along the z axis 
in the field of a helical undulator. In what follows we use the following as- 
sumptions: i) the electrons move along constrained helical trajectories in par- 
allel with the z axis; ii) the radius of the electron rotation in the undulator, 
r w — @w/k w , is much less than the transverse size of the electron beam. Next 
let us assume that electron beam density at the undulator entrance is simply 
n = n (r*j_)[l + a in cos ou(z/v z — t)}, where a in = const. In other words we con- 
sider the case in which there are no variations in amplitude and phase of the 
density modulation in the transverse plane. 

Under these assumptions the transverse current density may be written in the 
form 

j ± = -ev ± (z)n(f ± ,z/v z -t) = -ev ± n (f ± )[l + a- m cosu(z/v z - t)] , 

where we calibrated the time in such a way that current density has its max- 
imum at time t — 0, at point z — 0. Even through the measured quantities 
are real, it is generally more convenient to use complex representation. For 

— * 

this reason, starting with real j±, one defines the complex transverse current 
density: 

jx + ij y = -ec9 w n (r ± ) exp(- i k w z)[l + a in cosu(z/v z - t)] . (1) 

Transverse current have the angular frequency uo and two waves traveling in 
the same directions with variations exp \{ujz/v z — k w z— out) and exp — i(uz/v z + 
k w z—ut) will add to give a total current proportional to exp(— i k w z) cos u(z/v z — 
t). The factor expi(uz/v z — k w z — out) indicates a fast wave, while the factor 
exp — \{uz/v z + k w z — out) indicates a slow wave. The use of the word "fast" 
("slow") here implies a wave with a phase velocity faster (slower) than the 
beam velocity. 
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Now we should consider the electrodynamic problem. Using Maxwell's equa- 
tions, we can write the equation for the electric field: 

c 2 V x (V x E) = -d 2 E/dt 2 - A-nd]/dt . 

Then make use of the identity 

V x (V x E) = V(V • E) - V 2 E , 



where V • E can be found from the Poisson equation. Finally, we come to the 
inhomogeneous wave equation for E: 

c 2 V 2 E - d 2 E/dt 2 = Anc 2 Vp + ATrdj/dt . (2) 



This equation allows one to calculate the electric field E(r, t) for given charge 
and current sources, p{r,t) and j{r,t). Thus, equation (2) is the complete and 
correct formula for radiation. However we want to apply it to a still simpler 
circumstance in which the second term (or, the current term) in the right-hand 
side of (2) provides the main contribution to the value of the radiation field. 
Since in the paraxial approximation the radiation field has only transverse 
components, we are interested in the transverse component of (2). Thus we 
consider the wave equation 

c 2 V 2 E ± - d 2 E ± /dt 2 = 4ndj ± /dt , (3) 



which relates the transverse component of the electric field to the transverse 
component of current density. 

We wish to examine the case when the phase velocity of the current wave is 
close to the velocity of light. This requirement may be met under resonance 
condition 

uj/c = uj/v z - k w . (4) 



First we may point out that the statement of (4), the condition for the relation 
between u, k w and v z , is the condition for synchronism between the transverse 
electromagnetic wave and the fast transverse current wave with the propagat- 
ing constant u/v z — k w . With a current wave traveling with the same phase 
speed as the electromagnetic wave, we have the possibility of (space) resonance 
between electromagnetic wave and electrons. If this is the case cumulative in- 
teraction between bunched electron beam and transverse electromagnetic wave 
takes place. We are therefore justified in considering the contributions of all 
the waves except the synchronous one to be negligible. 
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It would be nice to find an explanation of resonance condition (4) which is 
elementary in the sense that we can see what is happening physically. Let 
us call the plane electromagnetic wave velocity c, which is of course, related 
to the phase constant by k = oo/c. Then the velocity of the beam relative 
to the electromagnetic wave is merely c — v z . With this relative velocity, the 
frequency of the electromagnetic waves as seen by the electrons is give by 
u e = (c — v z )k, where we have used the obvious relation between the phase 
constant k and the wavelength of the electromagnetic wave, k = 2n/\. Re- 
lation for oj c can be written as uj c = uj — kv z . Under the condition when (4) 
holds, we see by comparison with that equation that this would be merely 
the undulator frequency u c = k w v z . Accordingly, the condition of (4) that we 
have imposed is merely that the electron velocity is such that as the electron 
run past the electromagnetic waves, the apparent frequency they see is the 
undulator frequency. This is the basic effect involved in this interaction. If 
the apparent frequency is the undulator frequency, a cumulative effect on the 
electromagnetic wave results, since the waves are being driven at a frequency 
corresponding to natural rotation frequency of the electrons in the magnetic 
field. Another way of saying the same thing is that in the frame of reference 
moving with the electrons, the electromagnetic waves have a Doppler-shifted 
frequency equal to undulator frequency. 

Any state of transverse electromagnetic wave can always be written as a linear 
combination of the two base states (polarizations). By giving the amplitudes 
and phases of these base states we completely describe the electromagnetic 
wave state. It is usually best to start with the form which is physically clearest. 

— * 

We choose the Cartesian base states and seek the solution for Ej_ in the form 



Here and in what follows, complex amplitudes related to the field are written 
with a tilde. The description of the field given by (5) is quite general. However, 
the usefulness of the concept of carrier wave number is limited to the case 
where the amplitude is slowly varying function of z. 

To determine the form of E Xjy (z : f±) we substitute (1) and (5) into (3). We 
have 



E. 



'x,y — 



E x , y (z,r±) exp[iu(z/c - t)] + C.C. 



(5) 



exp[iu(z/c-t)] \ Vi + 



2iu d d 2 \ (E x 
c dz dz 2 J \E y 



) 



+ C.C. 




e9 v/ a in no(f±) smu(z/v z 



t). 



(6) 



Here 



is the Laplace operator in transverse coordinates. 
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Now we have an apparently simple pair of equations - and they are still exact, 
of course. The question is, how to solve these equations? First we simplify 
the equations obtained by noting that for a radiation field it is reasonable 
to assume that E Xjy (z,r±) are slowly varying functions of z (because of the 
radiation directivity) so that d 2 E X:y /dz 2 may be neglected. The correspond- 
ing requirement for the complex amplitude is | dE x ^ y /dz |<^C k | E XjV |. In 
other words, the radiation pulse must not change significantly while traveling 
through a distance comparable with the wavelength A = 2ir/k. This assump- 
tion is not a restriction. Such is the case in all practical cases of interest. 
Differential equations becomes 



exp[iuj(z/c-t)} \ Vi + 



2 i ijj d 



c dz 



E x 
,E„ 



+ C.C. 



u ( cos(/c w z)\ (.,..(. , s 

c \ — sm{K w ) I 



(7) 



The reader may well wonder why such a transformation is useful; therefore we 
digress temporarily to address this question. 

Such a transformation immediately modifies the hyperbolic wave equations to 
the parabolic equations. We should stress that the equations (7) are preferable 
for an analytical solution, since the mathematical techniques are always con- 
nected with more conventional parabolic equations. Although equations (7) 
cannot be solved in general, we will solve them for some special cases. These 
equations can be simplified further by noting that the complex amplitudes 
E x>y will not vary much with z, especially in comparison with the exponential 
terms exp(— i k w z). The slow wave of transverse current oscillates very rapidly 
about an average value of zero and, therefore, does not contribute very much 
to the rate of change of E xy . So we can make a reasonably good approximation 
by replacing these terms by their average value, namely, zero. We will leave 
them out, and take as our approximation: 

^(t) = -(i>*>-^)exp(- iC .) . ,8) 



Even the remaining terms, with exponents proportional to C — u/v z —u/c—k w 
will also vary rapidly unless C is near zero. Only then will the right-hand side 
vary slowly enough that any appreciable amount will accumulate when we 
integrate the equations with respect to z. The required conditions will be met 
if 

C < K 1 < k w z . 
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In other words, we use the resonance approximation here and assume that 
complex amplitudes E XiV are slowly varying in the longitudinal coordinate. By 
"slowly varying" we mean that | dE x>y /dz |<^C k w \ E x , y |. For this inequality 
to be satistied, the spatial variation of E x>y within an undulator period A w = 
2ir/k w has to be small. 

Equations (8) are simple enough and can be solved in any number of ways. 
One convenient way is the following. Taking the sum and the difference of the 
two we get 



^ + ( E x + iE y ) = 2-Ki^-e6 w a in n (r ± )exp(-iCz) , 



(9) 



vi + ^^](^-i4) = o. (10) 



These equations describe the general case of electromagnetic wave radiation 
by the bunched electron beam in the helical undulator. Equation (9) and (10) 
refer to the right- and left-helicity components of the wave respectively. The 
solutions for the right- and left-helicity waves are linearly independent 1 . It 
follows from (9) and (10) that only those waves are radiated that have the 
same helicity as the undulator field itself. 

Of course we could predict such a result. The electric field, E±, of the wave 
radiated in the helical undulator in resonance approximation is circularly po- 
larized and may be represent in the complex form: 

E x + iE y = E(z,r ± )exp[iu(z/c-t)] . (11) 



Finally, the equation for E can be written in the form 

(vi + ^J^) E = 2ni^ee w a hl n (r±)exp(-iCz) . (12) 



Equation (12) is an inhomogeneous parabolic equation. Its solution can be 
expressed in terms of a convolution of the free-space Green's function (impulse 

1 In the general case these components are not independent, but are connected by 
the boundary conditions on the vacuum chamber walls. It is usually assumed that 
the vacuum chamber walls are placed far enough from the electron beam, formally 
at infinity. Such an approximation well describes FEL amplifiers operating in the 
visible down to X-ray wavelength range. To describe FEL amplifier operating in 
the far-infrared wavelength range one should take into account the influence of the 
walls on the amplification process 
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response) 



G(z-z',r ± -r* ± ) 



1 



4ir(z — z'\ 



exp 



2c(z - z') 



(13) 



with the source term. When the right-hand side of (12) is equal to zero, it 
transforms to the well-known paraxial wave equation in optics. 

The radiation process displays resonance behaviour and the amplitude of elec- 
tric field depends stronly on the value of the detuning parameter C. With the 
approximation made in getting (12) the equation can be solved exactly, but 
the work is a little elaborate, so we won't do that until later when we take 
up the problem of output radiation characteristics. Now we will find an exact 
solution for the case of perfect resonance. When the parameters are tuned to 
perfect resonance, such that C = 0, the solution of the equation (12) has the 
form 



E(z,r ± ) 



ie9 w ua- m f dz' 



2c 



J - ^ , J dr^rio^exp 



2c(z - z') 



, (14) 



where (z, r±) and (z 1 , r^_) are the coordinates of the observation and the source 
point, respectively. 

Let us consider an axisymmetric electron beam with gradient profile of the 
current density. In this case we have — ev z no(r±) = —j S(r), where r is the 
radial coordinate of the cylindrical system (r, 0, z) and S(r) describes the 
transverse profile of the electron beam. To be specific, we write down all the 
expressions for the case of a Gaussian transverse distribution: 

-ev z n (f ± ) = -j S(r) = -— ^ exp 




where Jo is the total beam current. At this point we find it convenient to 
impose the following restriction: we focus only on the radiation see by an 
observer lying on the electron beam axis. When r = and the beam profile is 
Gaussian, we can write (14) in the form 

„ z . oo 

o o 



(r'f 




2a 2 


exp 



A 2 



2c(z - z') 



dr' . 
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Integrating first with respect to r', we have 



E (z) 



i 6 w (jja- m I f 



dz' 



(15) 



2c 2 J 



z — z' + i ka 2 ' 



o 



where k = u/c is the radiation wavenumber. The physical implications of 
this result are best understood by considering some limiting cases. It is con- 
venient to rewrite this expression in a dimensionless form. After appropriate 
normalization it is a function of one dimensionless parameter only: 



where z, = z/L w is the dimensionless coordinate along the undulator, L w is 
the total undulator length, (3 = ka 2 /L w is the diffraction parameter and 

E = E/E n = 2c 2 E/(6 w uja m I ) 

is the normalized field amplitude. The changes of scale performed during the 
normalization process, mean that we are measuring distance along the un- 
dulator and electron beam size as multiples of "natural" undulator radiation 
units. Let us study the asymptotic behaviour of the field amplitude at large 
values of the diffraction parameter (3. In this case (3 ^> z — z' and we have 
asymptotically: 



Now let us study the asymptote of a thin electron beam. In this case (3 — > 
and the function /(£, (3) can be evaluated in the following way. First we remark 
that the integral can be expressed as p + i q: 



The first integral should go from to z/(3, but is so far from z/(3 that 
the curve is all finished by that time , so we go instead to oo - it makes no 
difference and it is much easier to do the integral. The integral is an inverse 
tangent function. If we look it up in a book we see that it is equal to tt/2. The 
second integral can be expressed as logarithm function. Thus we have 




E — > z/(3 as (3 — > oo . 



(16) 




E Q ^n/2 + \\n(z/(3) 



as (3 -> . 



(17) 
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It is interesting to note that in the case of a thin electron beam, the radia- 
tion field, acting on the electrons is almost constant along the undulator axis, 
which is a little strange, because the modulated electron beam radiates elec- 
tromagnetic waves along the undulator. But that is the way it comes out - 
fortunately a rather simple formula 2 . 

Special attention is called to the fact that in the thin electron beam case, at 
/3 — > 0, amplitude E(z) is a complex function. One immediately recognizes 
the physical meaning of the complex E (z). Note that electric field (response) 
is given by fast wave of transverse current ("force") times a certain factor. 
This factor can either be written as p + i q, or as a certain magnitude p times 
exp(i5). If it is written as a certain amplitude p times exp(i<5), let us see what 
it means. This tells us that electric field is not oscillating in phase with the fast 
wave of transverse current, which has (at C — 0) the phase ip = uz/c — out, 
but is shifted by an extra amount S(z). Therefore S(z) represent the phase 
shift of the response. 

To get an intuitive picture on what happens with the radiation beam inside 
the near zone according to equation (14), let us first choose thin beam asymp- 
totic. This is an example in which diffraction effects play an important role. 
Simple physical consideration can lead directly to a crude approximation for 
the radiation beam cross-section. There is a complete analogy between the 
radiation effects of the bunched electron beam in the undulator and the radi- 
ation effects of a sequence of periodically spaced oscillators. The radiation of 
these oscillators always interferes coherently at zero angle with respect to the 
undulator axis. When all the oscillators are in phase there is a strong intensity 
in the direction 6 = 0. An interesting question is, where is the minimum? If 
we have a triangle with a small altitude r ~ z6 and a long base z, then the di- 
agonal s is longer than the base. The difference is A = s — z ~ r 2 /2z ~ z6 2 /2. 
When A is equal to one wavelength, we get a minimum. Now why do we get 
a minimum when A ~ A? Because the contributions of various oscillators are 
then uniformly distributed in phase from to 2tt. In the limit of small size of 
the electron beam interference will be constructive within an angle of about 
9 C ~ 1/Vkz. 

In the limit of large electron beam size, the spread of the angle distribution can 
be found from a two-dimensional Fourier transform of the field. The radiation 
field across the electron beam may be presented as a superposition of plane 
waves, all with the same wavenumber k = u/c. The value of k±/k gives the sine 
of the angle between the z axis and the direction of propagation of the plane 
wave. In the paraxial approximation k±/k = sin# ~ 9. We can expect that 



2 For FEL experts who happen to be reading this we should add that our logarith- 
mic term in (17) and the logarithmic growth rate asymptote for conventional FEL 
amplifier at small diffraction parameter (see [12]) are ultimately connected 
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apparent 




L w 

source length 

Fig. 4. Apparent photon source size in the case of a thin electron beam. The photon 
source extends over undulator length L w along the particle path 

the typical width of the angular spectrum should be of the order 9 C ~ (ka)^ 1 , 
a consequence of the reciprocal width relations of the Fourier transform pair 
Ak±a ~ I. 

The boundary between these two asymptotes is about ka ~ \/kL w or (another 
way to write it) a 2 ~ <rj if = L w jk. A rough estimate for the diffraction effects 
to be small is a 3> L w / (ka), which simply means that the diffraction expansion 
of the radiation at undulator length must be much less than the size of the 
beam. Another way to write this condition is ka 2 /L w — (5 ^> 1. The parameter 
(3 can be referred to as the electron beam Fresnel number. 

Our estimates show (see Fig. 4) that in the case of a thin beam the square of 
the radiation beam grows linearly with the undulator length 

S rad ~ {z6 c ) 2 ~ z/k for < 1 . 

In the case of a wide electron beam the most of the radiation overlaps with 
the electron beam 

S'rad — a 2 = const. for 1 <g P . 

Up to this point we only talked about electric field amplitude on the electron 
beam axis. Ultimately we want to consider the distribution of the electric 
field over the full electron beam cross-section. The full analysis is a difficult 
analytical problem. The calculations can be performed without great difficulty 
in two limiting cases, namely, the cases of diffraction parameter very large and 
very small comared with unity. To calculate equation (14) we note that the 
behaviour of Green's function (13) for ka 2 /L w ^> 1 approaches the behaviour 
of the delta function. The source function n (r^_) does not vary very much 
across the region | f i — | 2 ~ L w /k in the case of a wide electron beam: 
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therefore we can replace it by a constant. In other words, we simply take 
n o{^±) outside the integral sign and call it n (f±). In this case the integral 
over source coordinates in (14) is calculated analytically 

E{z, f±) = ire9 w za in n (r±) = ^^2° exp (~^2 ) as ^ ^ 1 ' ( 18 ) 



We emphasize the following features of this result. At exact resonance (C = 0) 
amplitude E is a real function. This tells us that electric field is oscillating in 
phase with the fast wave of transverse current. 

Let us turn next to the second case where (3 <C 1. It should be noted that 
dealing with klystron amplifier this situation is more preferable than the pre- 
vious one. In thin electron beam case there is only a small variation in electric 
field amplitude over electron beam cross section. This feature follows natu- 
rally from the fact that the radiation expands out and electron beam becomes 
thin with respect to the radiation beam. Using (14), we find that the field 
amplitude inside the thin electron beam is given by the formula 



E(z,r±) = 



2c 2 



71 , / CZ 

- + iln — 

2 Vcucr 



as < 1 . (19) 



The nonzero imaginary term is a manifestation of the fact that electric field 
is not oscillating in phase with the fast wave of transverse current. 

In (14) we have the general solution of the problem of finding the field produced 
by modulated current, and this is the first part of the problem of the operation 
of the high gain klystron amplifier. The next question which we must take up 
is the converse one: Given an electric field of the type we have just considered, 
what is the electronic motion in such a field? 



2.2 Self -interaction within a modulated electron beam moving in an undulator 



Since, the wavelength of the radiation pulse is much less than either the trans- 
verse electron beam size or the pulse length it can be accurately represented 
by a wave traveling at u/k ~ c. Further, since the electron also travels at 
v z ~ c, the electron and the radiation remain essentially overlapped during 
the transit through the undulator. Although the radiation fields alone do not 
affect the electron trajectory, the interaction of the transverse radiation elec- 
tric field with the transverse oscillations of the electron induced by the static 
undulator magnetic field causes the electron to either gain or lose energy. 

As simple example of the interaction between electrons and electromagnetic 



20 



wave in an undulator we shall consider a circularly polarized plane electromag- 
netic wave propagating parallel to the electron beam. The field of the elec- 
tromagnetic wave has only a transverse component, so the energy exchange 
between the electron and the electromagnetic wave is due to the transverse 
component of the electron velocity. The rate of electron energy change is 

d£ 2 d7 - £? 
= m c c — = — ev j_ • ti\ , 

dt dt 

— * 

where E± is the vector of the electric field of the wave: 
E±_ = E {e x cos [uj{z/c — t)\ + e y sin [u(z/c — t)]} . 

Remembering that d z = v z d t we find 

— = [V X E X + VyEy) 

dz v z 

~ —e9 w E {cos(k w z) cos [u(z/c — t)] — sin(k w z) sin [u(z/c — t)]} 
= —e9 w E cos [k w z + uo(z/c — t)\ 

= -e6 w Ecos*p. (20) 

The phase tp has a simple physical interpretation and is equal to the angle 
between the transverse velocity of the particle, v±, and the vector of the elec- 
tric field. For effective energy exchange between the electron and the wave, the 
scalar product (ev±-E±) should be kept nearly constant along the whole undu- 
lator length, i.e. a synchronism should be provided. This resonance condition 
may be written as 

d^ = k w dz + (u/c) dz - codt = . 

Remembering that dz = v z dt, we have k w + u/c — lj/v z = 0. Or, since 

A = 2-kc/uj, 

K/Vz = \/(c-V z ) . 

Thus, we see that synchronism takes place when the wave advances the elec- 
tron beam by one wavelength at one undulator period. This resonance con- 
dition is exactly equation (4), so we have come full circle. Electromagnetic 
wave excitation and production of the electron beam energy modulation are 
ultimately connected. 

From this simple analysis it is apparent that the electron's relative position 
within a radiation wavelength will determine whether it consistently gains or 
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loses energy as it travels through the undulator magnet. A convenient variable 
for the description of the interaction between electrons and electromagnetic 
wave in an undulator is the phase ip. The relevant value of the phase ip is that 
at the location of the particle. Hence, the total derivative of ip is given by 

d ^ _ £, _ dip dip dt _ ^ + ^ oj 
dz dz dt dz w c v z 

Thus, the phase of the particle changes when the resonance condition is not 
satisfied exactly. 

Since discrete electrons are considered, their entrance into the interaction re- 
gion is conveniently described in terms of their entrance phase positions rel- 
ative to one cycle of the modulating electromagnetic wave at the input, i.e. 
z — 0. An entrance phase variable is defined by ipo = —cuto, where to is the 
entrance time. Let the initial electron energy be So- Thus, at a given displace- 
ment plane the electron will have a kinetic energy deviation 

z 

£(ip , z) - £ = ~e9 w J dzE cos (ip + Cz) , 

o 

if ipo is the arrival phase of this particular electron at the undulator entrance. 

We should remark that our analysis of the interaction between electrons and 
radiation in an undulator gives a result that is somewhat simpler than we 
would actually find in nature. All of our calculations have been made for the 
plane electromagnetic wave with phase velocity equal to the velocity of light. 
The actual radiation wave in an undulator will experience an amplitude change 
and a phase shift due to the radiation process. The electric field of the wave 
radiated in the helical undulator can be represented by equation (11). We can 
also write this expression as 

E(z : f±)expiuj(z/c-t) =\ E(z,r±) \ expi[uj(z/c — t) +ip e (z,f±)] , 

which says that radiation wave is not ordinary transverse plane wave. Our 
final result for the kinetic energy deviation is therefore 

z 

£ — So = —e9 w J dz \ E(z, f±) | cos [ip + Cz + ip c (z, f±)] , (21) 
o 

This expression deserves some remarks. In dealing with start-up from the elec- 
tron beam density modulation the appropriate "frame of reference" is moving 
with the fast wave rather than at the velocity of light. The angle between the 
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transverse velocity of the particle and the vector of the electric field is given 
by 



where Cz + ip e (z : f±) represent the phase shift of the radiation wave ("re- 
sponse") relative to the fast wave ("force"). What do we use for the entrance 
phase in our formula (21)? We use the entrance phase position relative to one 
cycle of the fast wave at the input. For simplicity of analysis, we consider the 
cases in which there are no variations in phase of the density modulation in 
the transverse plane at the undulator entrance. In those cases an entrance 
phase variable is defined by ip = —ut , where t is the entrance time. 

We have what we need to know - the electron beam energy modulation at the 
undulator exit (21). We are ready to find energy modulation, because we have 
already worked out what field is produced by a bunched electron beam inside 
the near zone (14). These equations are not well suited for quick calculation of 
the energy modulation at a particular diffraction parameter. We may, however, 
express (14) in much simpler form for very small and very large diffraction 
parameters, making use of limiting expressions (18) and (19). Let us see what 
happens if the diffraction parameter (3 is large. At exact resonance (C = 0) 
using (18) the energy modulation amplitude achieved at a given displacement 
plane can be written as: 

S£ <*Ji z* Jo / r 2 \ 



where I a = m c c 3 /e ~ 17kA is the Alfven current. 

Now let us analyze this expression and study some of its consequences. To 
develop a useful mathematical formalism for the description of interaction 
between electrons and fields in an undulator we must choose the appropriate 
coordinate system to maximize physical clarity. We chose at the beginning of 
this section to use the stationary coordinate system z. For example, the beam 
density at point z, at time t is 

n(z/v z — t) = n + a in n cosu(z/v z — t). 



A convenient way of description the particle beam dynamics is to use a co- 
ordinate system z' that moves at a velocity v z with respect to the stationary 
coordinate system z. When a particular electron enters the undulator, its ini- 
tial position is z' in the moving coordinate system. The electron would remain 
stationary at z' in the moving coordinate system; in the stationary coordi- 
nate system, at a time t, it would have moved a distance from the entrance 
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z = v z (t—to), where to is the entrance time. Let us write an expression for beam 
density modulation in terms of the moving coordinate system introducing the 
wave number k' defined by k' — uj/v z — 2ir/\', where A' is a wavelength of the 
modulation frequency, i.e. the distance through which a particle travels in a 
cycle of modulation frequency. We may thus write n(z' ) = n + am^o cos (k'z' ). 
Note that k'z' = —ut = tp and this expression may be rewritten in the terms 
of arrival phase 

n(ifj ) =n + a in n cost/> , 



where we calibrated time in such a way that when a particular electron enters 
the undulator at time t — 0, its initial position is z' = in the moving 
coordinate system. Just as we expected, in the case of wide electron beam, the 
beam energy modulation is oscillating in phase with beam density modulation. 

In the case of thin electron beam, things are quite different. For a small diffrac- 
tion parameter < 1 we may apply an asymptotic approximation for the field 
amplitude (19) and get 



58_ 



\ujz Iq 



2c >yI A 



71 i 
2 cosVo 



In 



cz 



sin ip 



(23) 



Thus, we see that, in the case of thin electron beam, the beam energy mod- 
ulation is not oscillating in phase with beam density modulation. Equations 
(22) and (23) give us an idea of what we should expect. Generally we can try 
to calculate the energy modulation precisely by using the (14) and (21). That 
is the end of our calculations of the electron beam energy modulation, but 
there is one physically interesting thing to check, and that is the conservation 
of energy. 



2.3 Power balance 



As any oscillating charge radiates energy, so must a modulated electron beam 
moving along an undulator radiate energy. If the system radiates energy, then 
in order to preserve conservation of energy we must find that electron beam en- 
ergy is being lost. An interesting question is, what forces are electrons working 
against? That is interesting question which can be completely and satisfacto- 
rily answered for modulated electron beam in an undulator. What happens in 
this: in bunched beam, the fields produced by the moving charges in one part 
of the undulator react on the moving charges in another part of the undulator. 
We can calculate these forces and find out how much work they do, and so find 
the right rule for the radiation force. We can make this calculation because at 
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short distance we know the electric field. Above we calculated the radiation 
field inside the near zone (14). 

In our system energy is transported both in the form of electromagnetic waves 
and in kinetic energy. The well-known Poynting vector represents the electro- 
magnetic power flow; the kinetic-power flow is merely the number of electrons 
crossing per unit area per unit time multiplied by the kinetic energy per elec- 
tron. Consider first the electromagnetic power. In the paraxial approximation 
the diffraction angles are small, the vectors of the electric and magnetic field 
are equal in absolute value and are perpendicular to each other. Thus, the 
expression for the radiation power, W, can be written in the form: 



where (• • •) denotes averaging over a cycle of oscillation of the carrier wave. If 
we consider a system with fields and bunched electron beam in an undulator, 
the energy stored in any volume fluctuates sinusoidally with time. But on the 
average there is no increase or decrease in the energy stored in any portion 
of the volume, so that the below conservation theorem holds if averaged over 
time, although the integrals might not be zero instantaneously because of the 
time variations in the energy storage. 

Since the radiation field inside the near zone has the form of equation (14), 
than W is given by: 



eV ! gag 1 f _ \ f d z' f _ 
= J d ^\J—J dr ^ oM 6XP 

\ f dz" f ^ _,, \ iu\r±-?i | 2 



The product of integrals over z' and z" can be represent as 

2 2 2 2' 

J <$>{z')&z' j $VW = J <t>(z')dz' J <£>*(z")dz" + C.C. 



The integral over transverse coordinate fj_ is equal to 
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As a result, expression (24) can be written in the form: 



W 



x exp 



j dz' j^j J dr\J dr> (rl)n (rl) 



I —i —if 1 9 

1 uj I rj_ — r_[ | 
2c(z' - z") 



+ C.C. 



(26) 



Latter expression can be written in terms of the complex amplitude of the 
radiated field E(z, f±). Using (14) and (26) we obtain the expression for W: 

w= ie^p L j dz > J d 7 ±no{ 7 ±) E(z>y ± ) + C.C. (27) 




The radiation power, W, must be equal to the difference of the electron beam 
power at the exit and the entrance of the undulator. The rate of the energy 
change for a single electron is given by 

d£/ dz = -e0 w I E(z,f ± ) | cos[^ + A(z,f ± )} . 



To obtain the mean power loss by the electron beam, we should multiply this 
value by the particle flux density v z a in n cos ipo ~ ca; n n cos ipo average over 
time, and integrate over the beam cross-section and the undulator length. 
Finally, we obtain the following result: 



AW C = ^r/ d ^o J dz' J dr*_L {ca hl n (r' A _) cos ip } 


x {-e6 w I E(z', r ± ) | cos^o + M^', ?±)}} 



26 



= _ie^ j dz >J djf±no ^ ±) E^y ± ) + c.C. (28) 
o 

Comparing this expression with (27), we see that the radiation power and the 
change in the electron beam power have equal absolute values and are opposite 
in sign, i.e. AW C + W = 0. So, power balance takes place. 



3 The region of applicability 

Our theory of klystron amplifier is based on the assumption that beam den- 
sity modulation does not appreciably change as the beam propagates through 
the undulator. When the resonance condition takes place, the electrons with 
different arrival phases acquire different values of the energy increments (pos- 
itive or negative), which results in the modulation of the longitudinal velocity 
of the electrons v z with the radiation frequency uo. Since this velocity mod- 
ulation is transformed into density modulation of the electron beam during 
the undulator pass, an additional radiation field exists because of variation 
in amplitude density modulation. Instead, we assume that the amplitude of 
the electron beam density modulation has the same value at all points in the 
undulator. This approximation means that only the contributions to the radi- 
ation field arising from the initial density modulation are taken into account, 
and not those arising from the induced bunching. 

3.1 Induced bunching 

It is interesting to estimate the amount of bunching produced during the 
undulator pass. The velocity v z of the electron at given displacement plane z 
is given by 

v z (£) = v + Sv z (z,r±) cos[^o + S(z,f±)] 

where vq = v z {£q) is the injected velocity, and 5v z is a small amplitude of high 
frequency oscillation. From (23) we can find Sv z remembering that v z ~ c and 
(dv z / d£) \s=£ ^ c/{~f 2 z £ ). We have 

5v z (z) ~ a in 9lcuzl (tzVa) _1 • (29) 

We now examine the behaviour of these electrons as they travel a certain 
distance in an undulator. Since successive electrons are not traveling with 
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the same velocity, bunching takes place. To calculate a numerical value for 
the current, we must compute the actual arrival time of the various electrons 
at the undulator exit. The appropriate expression is formulated by means of 
the law of conservation of charge. Quantitatively, we can say that during an 
infinitesimal interval of time cUi, a total charge dq = hdt\ passes through 
plane z±. We follow the same particles (i.e. the same charge) through the plane 
Z2 at some later time t 2 . These particles will pass through the z 2 plane during 
the an infinitesimal time interval d^2 which may be shorter or longer than 
cUi- In either case, the instantaneous current I 2 (t 2 ) due to the charge passing 
through z 2 plane in the interval d^2 is obtained from dq = h dt 2 . Obviously, 
the two charges are the same since we are following the same set of particles. 
Hence due to this charge, there is a current through Z2 plane I2 — h dti/ dt 2 . 
If an electron has velocity v z when it leaves point z\ and maintains this velocity, 
then its transit time to point z 2 is given by the expression 



t 2 = * 1 + 



Z2 ~ Zl 



Vq + Sv z cos uti 



tl + 



Z2 - Z X 



1 coswti 

^0 



(30) 



This equation provides us with the relation between the arrival time i 2 at point 
z 2 for the particular group of charges which left point Z\ at the departure time 
t\. It is possible to calculate current at point z 2 for small values of Sv z /v by 
differentiating the last equation and substituting in equation I 2 = hdti/ dt 2 . 
This procedure leads to 

h ^ — — —. — ^ h (1 - asin^ti) , 

1 + a smwti 

where a = uj(z 2 — Z\)8v z jv\. This equation is obviously still not in a useful 
form, since it also only gives the current at point Z2 in terms of departure 
times of the charges from point Z\. However, by using (30) again (to eliminate 
ti) and using the fact that a C 1, one can write approximately 

I 2 ~ h {1 - ccsin [ut 2 — k(z 2 — ^i)]} , 



where the terms that have been neglected in t\ would have contributed only 
terms of the order of a 2 to the current Thus the induced bunching that we 
want at the point z 2 due to the velocity modulation at the point z\ is equal 
to a = ijj{z,2 — zi)Sv z /v 2 . Now we find the total induced bunching at the point 
z 2 by finding the induced bunching there from each point Z\ < z 2 , and then 
adding the contributions from all the points. To calculate this sum we need 
to use (29). We have, then, 

Sa = -2 / (z 2 - zi) —^dzi ~ a m uj 2 9lzll (c 2 ^I A ) ■ (31) 
v J dzi v ' 
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Thus, the requirement for the induced bunching to be small can be written as 
5a <C a in . The result of a more careful analysis shows that the last condition 
can be written as 

c 2 ^/o(6c 2 7 2 V A ) _1 «l. (32) 



3.2 Energy spread 



During the passage through an undulator the electron density modulation can 
be suppressed by the longitudinal velocity spread in the electron beam. For 
effective operation of klystron amplifier the value of suppression factor should 
be close to unity. If there is an (local) electron energy spread into the electron 
bunch, AS, there will be a corresponding longitudinal velocity spread given by 
cAS I (7 2 £o) • To get a rough idea of the spread of electron density modulation, 
the position of the particles within the electron beam at the undulator exit 
has a spread which is equal to 

Az' ~ L w Av z /v ~ L W AS (7^0) ~* , 



We know that uncertainty in the phase of the particles is about Aip ~ 
2itAz'/\'. Therefore, a rough estimate for the microbunching spread to be 
small is 

uL w AS (c-ilS^y 1 < 1 . 



It is clear that we have made several rough approximations. Besides various 
factors of two we have left out, we have used AS, where we should have 
used variance {(AS) 2 ). All of these refinements can be made; the result of a 
more careful analysis shows that for the Gaussian distribution function the 



suppression factor is equal to exp — u; 2 L^((A£) )/ (2c 2 7^£q) 



So a better 



answer is 



<(A£)' 

2c 2 7 



S 2 



< 1 



Remembering that A w = 27 2 A and v ~ c, we have 
(2vriV w ) 2 ((A£) 2 >/£ 2 «i , 



(33) 



where N w is the number of undulator periods. This condition is not a re- 
striction. Let us present a specific numerical example. For unbunched electron 
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keV. If the nominal 



energy is £ = 500 MeV, the last condition can be written as N w <C 10 4 . Such 
is the case in all practical cases of interest. 



3.3 Transverse emittance 

For a klystron amplifier it is of great interest to maximize the gain per cascade 
pass which is proportional to the amplitude of electric field inside the electron 
beam. It is usually desired to optimize the amplifier to achieve the highest 
possible gain coefficient for a given total current Jo- Reducing the particle beam 
cross-section by diminishing the betatron function reduces also the size of the 
radiation beam and increases, according to (18), the amplitude of the electric 
field. This process of reducing the beam cross section is, however, effective 
only up to some point. Further reduction of the particle beam size would 
have practically no effect on the radiation beam size because of diffraction 
effects (19). It is therefore prudent not to push the particle beam size to 
values much less than the difraction limited radiation beam size. From the 
preceding discussion we may want to optimize the beam geometry as follows. 
For maximum amplification, the transverse size of the electron beam has to 
be matched to the diffraction limited radiation beam size 



The wavelength and undulator length dictates the choice of the optimal trans- 
verse size of the electron beam. 

One might think that all we have to do is to get electric field amplitude to 
a maximum - we can always reduce electron beam cross section and we can 
always increase gain. So it is not impossible to build an electron focusing sys- 
tem that has a small betatron function. In fact, one of the main problems in 
klystron amplifier operation is preventing the spread of microbunching due to 
angle spread into the electron beam. Some electrons traverse the undulator 
not along or parallel to the z— axis, but at a small angle 6. If there is an rms 
angular divergence a' within the electron bunch there will be a correspond- 
ing longitudinal velocity spread given by Av z /v z ~ (a') 2 /2. We assume that 
the amplitude of beam modulation does not appreciably change as the beam 
propagates through the undulator. To avoid an decrease in the amplitude of 
the beam modulation, the longitudinal velocity spread must be small enough 
L w Av z /v z <C X/(2tt). This gives the condition 




(a') 2 < A/(2vrL w ) . 
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The last two conditions can be written as 
e/vr = aa' < A/(2tt) , 



where e is the transverse electron beam emittance. The last condition can be 
seen as "phase matching" of the electrons and radiation. To produce maximal 
gain the particle beam emittance must be less than the diffraction limited 
radiation emittance. Obviously, this condition is easier to achieve for long 
wavelengths. For optimal gain for a 30 nm radiation, for example, the electron 
beam emittance must be smaller than about 10~ 6 7r rad. — cm. This condition 
may be easily satisfied in practice. 



4 Bunching by energy modulation 

In two-undulator klystron, energy modulation occurs in the first undulator and 
the conversion to density modulation in the dispersion section. The dispersion 
section is designed to introduce the energy dependence of a particle's path 
length, Az = R 5 q5£/£o. Several designs are possible, but the simplicity of a 
four-dipole magnet chicane is attractive because it adds no net beamline bend 
angle or offset and allows simple tuning of the momentum compaction, i? 56 , 
with a single power supply. The trajectory of the electron beam in the chicane 
has the shape of an isosceles triangle with base length L. The angle adjacent 
to the base, 0b, is considered to be small. For ultra-relativistic electrons and 
small bend angles, the net _R 56 of the chicane is given by (see Fig. 5) 

-^56 = L9ft . 

Optimal parameters of the dispersion section can be calculated in the following 
way. The phase space distribution of the particles in the first undulator is 
described in terms of the distribution function f(P,ip ) written in "energy- 
phase" variables P = £ — £ and ip = 2nz'/\' = uo(z/v — t), where £ is the 
nominal energy of the particle and uj is the angular frequency. Before entering 
the first undulator, the electron energy distribution is assumed to be Gaussian: 

/o(P) = 75^) exp ("2«&)) ' 

The present study assumes the density modulation at the end of first undulator 
to be very small and there is an energy modulation P cos ip only. Then the 
distribution function at the entrance to the dispersion section is 

/o(P + P sinVo) • 
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Fig. 5. Schematic of dispersion section 



After passing through the dispersion section with dispersion strength d VW dP, 
the electrons of phase ip and energy deviation P will come to a new phase 
^0 + P d^Po/ dP- Hence, the distribution function becomes 



f(P,ip )=fo P + PosinUo 




The dispersion strength parameter and compaction factor are connected by 
the relation 

d^o _ 27rdy _ 27ri? 56 
dP ~ \> d£~ A' £0 ' 



The integration of the phase space distribution over energy provides the beam 
density distribution, and the Fourier expansion of this function gives the har- 
monic components of the density modulation converted from the energy mod- 
ulation [13]. At the dispersion section exit, we may express current I in the 
form 



~ 00 
I = I j /(Pi„)dP = /o + 2/„^exp 

—00 

d^o 



71=1 



x J n nP n 



dP 



cos(nifj Q ) . 



(34) 



We find a set of harmonic waves, of which the fundamental term, with angu- 
lar frequency u, is the one of importance in an amplifier. This fundamental 
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involves the phase variation cos^, but it is multiplied by the amplitude term 



a out = 2Ji I -n)-j-p I exp 



For small input signal we may assume that the argument of the Bessel func- 
tion is small. The function Ji(X) approaches X/2 for small X; thus the mi- 
crobunching approaches 



flout 



Po^exp 



We see that microbunching depends greatly on the choice of the dispersion 
section strength. An attempt to increase of the amplitude of the fundamental 
harmonic, by increasing the strength of dispersion section, is countered by 
a decrease of the exponential factor. The microbunching a out has clearly a 
maximum 

(°out) m = = (35) 

/2.72((A£) 2 > 



and the optimum strength of the dispersion section is 

/djM i 

UWma* x/<(A£) 2 > ' 



(36) 



Consider the situation at ^((A^) 2 ) = 5 keV, £ = 500 MeV, A = 30 nm. 
Appropriate substitution in (36) show that the optimal compaction factor is 
equal to R 56 = LQ\ = 500 /im. 



5 The properties of a klystron amplifier 



We can now piece together the information we have obtained in the preceding 
sections, and discuss the behaviour of the klystron as an amplifier of the 
electron beam intensity modulation. The present section is concerned with 
cascade amplifiers designed for conditions approximating maximum gain per 
cascade pass. 
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5.1 Klystron amplifier gain 



In the section 2 we learned that when a modulated electron beam passes 
through an undulator, the radiation field modulates the energy of the elec- 
trons. In order to calculate the energy modulation we can use (14) and (21) 
and note that energy modulation at the undulator exit generally can be writ- 
ten as pcosips + <?sin-?/>o, or as pcos(-^o + 5). When the gain per cascade pass 
is high enough (a ou t din), the cascades can be considered independent. In 
this case we do not see the effect of phase shift 5, but see only total amplitude 
of modulation equal to p = \/p' 2 + q 2 . 

We are going to apply (14) to our analysis of energy modulation taking ad- 
vantage of similarity techniques we used in section 2.1. The beam, of peak 
current Jo, is modulated in energy by the input undulator, and acquires an 
energy modulation amplitude on the electron beam axis, of amount 




where 




0; n 



2c 7 / A ' 



and universal function F((3) should be calculated numerically by using (14). 
The physical implication of this result are best understood by considering 
some limiting cases. Equations (22) and (23) give us an idea what we should 
expect in the case of wide and thin electron beam. We have asymptotically: 

F{(3) -> 1/(2/5) as (3 -> oo , 

F{(3) -> ^tt 2 /4 + Ml/P) - I] 2 as (5 -> . 



A natural and interesting choice is to calculate the gain in thin beam case 
which we can identify with maximal gain. Note that the dependence of the 
factor F((3) on the exact size of the bunch in the thin beam asymptote is 
rather weak. Thus F(@) ~ 1 is a reasonable approximation. Finally, the energy 
modulation amplitude can be estimated simply as: 

58 ( 58\ a-mOtuLv, I n , _ , _, 

for (3 < 1 . (37) 



8 \8L 2c 7 ^ 



The energy modulated beam then enters the dispersion section, in which den- 
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sity modulation is developed. From equations (35) and (37) the over-all gain 
per cascade pass at resonance may be readily be obtained if optimum disper- 
sion section strength is assumed; thus 

rt a out ^W L W Iq . , 

G = ~ — , - — — for p < 1 , (38) 



where the following notation has been introduced: A7 = ^ {(AS) 2 ) / m 2 c 4 . 
Putting 6 W = K/i, L w = N W X W , uj = 2ir/\ = 4n'y 2 c/X w and 7^ = 7 2 /(l + ^ 2 ) 
in (38) gives 

K 2 N T 

G °^(TT7^ for ' 3<1 - (39) 



The equation (39) tells us the maximal gain of cascade as a function of number 
of undulator periods N w , undulator parameter K, peak current I and rms 
local energy spread A 7 . 

For conventional XFELs the following problem exists. The required very small 
transverse beam emittance can be obtained only in a photoinjector. The bunch 
length must be shortened. The distribution of particles in phase space is given 
either by the injector characteristics and the injection process. Longitudinal 
phase space can be exchanged by special application of magnetic and RF 
fields. This is done in a specially designed beam transport line consisting of a 
nonisochronous transport line (magnetic chicane) and an accelerating section 
installed at the beginning of the bunch compressor. Liouville's theorem is not 
violated because the energy spread in the beam has been increased through 
the phase dependent acceleration in the bunch compression system. 

The cascade gain with continually optimized strength of the dispersion section 
would be directly proportional to Iq/A^. We have a rather surprising result. 
We know that ratio Iq/A'j depends on the parameters of photoinjector but 
not on the bunch compression. We expect this ratio to depend only on the 
longitudinal emittance. So the simple result says that the klystron gain is 
independent of absolute value of peak current. The change in peak current 
with change of bunch length is compensated by the larger energy spread. Since 
/0/A7 is independent of bunch length, the gain would then be independent of 
bunch compression. 

The formula for the gain which we derived (39) refer to the case of the helical 
undulator. For somewhat wider generality, although we are still making some 
special assumption about undulator magnetic structure, we shall calculate the 
characteristics of the klystron amplifier with a planar undulator. The magnetic 

— * 

field on the axis of the planar undulator is given by H = e x H w cos (k w z). The 
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explicit expression for the electron velocity in the field of the planar undulator 
has the form: v y = —e y c9 w sm(k w z), where 9 W = K/j = \ w eH w / ' (2imi e c 2 ^) . It 
is not hard to go through the derivation of electron beam energy modulation 
again. If we do that, and calculate the gain the same way, we get 

Go = ^ ~ 9 l^h \ for p < i , (40) 
a in 4 v / 27T2c A 7 / a h y ' 



where 

Au = [J (Q) - Ji(Q)} , 

J n (Q) is a Bessel function of nth order, 
Q = ^> /(8A; w 7 2 ) = K 2 /(A + 2K 2 ) . 



When we simplified the expression for Q, we used the resonance condition for 
the planar undulator ujq = 2'~f 2 k w /[c(l + K 2 /2)]. It is convenient to rewrite the 
expression (40) in the form: 

G »-^(#^)^ for < 41 » 



The significance of the proposed scheme cannot be fully appreciated until we 
determine typical values of the gain per cascade pass that can be expected in 
practice. Let us present a specific numerical example for the case of a klystron 
amplifier with a planar undulator. With the numerical values A w = 3 cm, K = 
1.42, 7 = 10 3 , the resonance value of wavelength is A = 30 nm. If the number 
of the undulator period is iV w = 100, normalized transverse emittance e n = 
2n urn, and betatron function is equal to the undulator length, the diffraction 
parameter is about (5 ~ 0.4. Remembering the results of section 2, we come to 
the conclusion that we can treat this situation as an klystron amplifier with 
thin electron beam. It is shown that the residual energy spread in the TTF 
FEL injector is on the order of a few keV only, even at a peak current of 100 
A. If I = 100 A, and \J ((AS) 2 ) = 5 keV, appropriate substitution in (41) 
shows that the intensity gain per cascade pass is about Gq ~ 10 4 . 

It has been seen above that there are definite upper limits to the gain of a 
two undulator klystron with a given current. It is obvious that a gain much 
higher than that allowed by these limits may be obtained by using more than 
one stage of amplification. In the three-undulator type of klystron, called the 
cascade-amplifier klystron, the additional undulator (and dispersion section), 
which lies between the input and output undulators, has no need of radiation 
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phase matching. Second cascade is excited by the bunched beam that emerges 
from the first cascade, and it produces further bunching of the beam. Cascade 
bunching in a high-gain cascade amplifier results in current components just 
like those of simple bunching; the equivalent intensity gain, which determines 
the cascade bunching, is given by G 2 C<1 = Gq™, where m is the number of 
cascades. In practice, with a two-cascade klystron amplifier, an intensity gain 
in excess of 80 dB may be obtained at wavelengths around 30 nm. This simple 
description is valid when beam modulation amplitude at the exit of the last 
cascade is much smaller than unity. What happens when this condition is not 
met is discussed later. 

Fluctuations of the electron beam current serve as an input signal to the 
XFEL. The initial modulation of the electron beam is defined by the shot 
nose and has a white spectrum. When the electron beam enters the first un- 
dulator, the presence of the beam modulation at frequencies close to resonance 
frequency initiates the process of radiation. The study has shown (see section 
5) that the mean square value of input is about (af n ) ~ 3 • 10~ 7 at 1$ = 100 A, 
N w = 100 and A = 30 nm. It is seen that saturation in the cascade klystron 
should occur at an exit of the second cascade in this situation. Thus, we con- 
clude that a klystron amplifier which consists of a succession of 2-3 cascades 
can operate as a soft X-ray SASE FEL. 

A disadvantage of using conventional FEL amplifier at small peak current is 
the enormously long undulator that is required. Let us calculate the param- 
eters of the conventional SASE FEL for the value I = 100 A of the peak 
current. Suppose yj '{(AS) 2 ) = 5 keV, e n = 2nfim, 7 = 1000, A w = 3 cm, 
K = 1.42, focusing beta function is equal to 3 m. Numerical solution of the 
corresponding 3-D eigenvalue equation (see [12]) demonstrates that the field 
gain length should be about l g ~ 3.7 m. Well-known that saturation in a 
SASE FEL with uniform undulator occurs after approximately 10 exponential 
field gain lengths. Remembering the total undulator length of the two-cascade 
klystron amplifier, we come to the conclusion that a significant impruvement 
in this parameter for klystron amplifier can be achieved. 



5.2 Klystron amplifier bandwidth 

All of the foregoing discussion of klystron amplifier gain has been concerned 
solely with the gain at resonance - that is u = u . Nothing has been said 
about amplification bandwidth. Now, we would like to find out how the gain 
varies in the circumstance that the seed signal frequency uj is nearly, but not 
exactly, equal to 47T7^c/A w . It is not hard to go through the derivation of 
electron beam energy modulation again. If we take C ^ 0, the solution of the 
wave equation (12) has the form 
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E(z,r ± ) 
x J dr'_ L n (r' 1 _) exp 



ie6 w ua in f dz' 



f—^-iC^) 

7 z — z 



2c J z — z 
o 



I -» -f |2" 

io; | rj_ — r ± | 
2c(* - 



(42) 



Using (21) the energy modulation achieved at nonzero detuning parameter 
can be written as 



58 = -iexp(i^ ) e2 ^ ain / d^exp(iCz') J exp(- i CV') 



dz" 



x / dr ± n (r ± ) exp 



o 

i a; | f j_ — fj_ 
2c(z' - z") 



+ C.C. 



(43) 



Now we shall calculate the energy modulation on the electron beam axis. 
When r = and beam profile is Gaussian, we have 



S£_ _ ifo aiJluj_Jo_ 



p iC(z'—z") 



Ac 'yip 



o o 



Let us see what happens if the diffraction parameter f3 is large. In this case 
we have asymptotically: 



Z Z 



£ Aa 2 jI A J J 



o o 



The integrals over z' and z" in the latter equation are calculated analytically 



z z 

J dz' j d/v^'-^ + c.c. 



o o 

sin 2 (Cz/2) , (Cz/2) - sm(Cz/2) cos(Cz/2) . , 

7 cos^o + - — ,„J„ — sm^o ■ 



(C/2) s 



(C/2) s 



(44) 



Taking into account the definition of the detuning parameter we find that C is 
connected by a simple relation with the frequency deviation: uj — u = Au = 
2^ 2 C. If we follow through the simple algebra we find that 
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Fig. 6. Frequency response curve 
6S_ ajJlLj Jp { sm 2 [nN w Au}/tJo} 
S 4 <7 2 7 Ja\ [ttN w Au/uj } 2 COS ^° 

[kN w Auj/ujo] — sin[7riV w Au;/u;o] cos[7riV w Acj/c<;o] . 
+ [nN w Au/u ] 2 81 



The energy modulation at the undulator exit can be written as pcostpo + 

/a 2 n is proportional to (p 2 + q 2 ). It is 



"out 



gsin^o- The intensity gain G 
interesting to plot this gain as a function of input signal frequency in order 
to see how sensitive it is to frequencies near the resonant frequency ojq. We 
show such a plot in Fig. 6, where the following notation has been introduced: 

x = ttN w Auj/ujo, 



G(x) 
G(0) 



N sin 4 (x) \x — sin(x) cosfa;)] 2 
j(x) = h 



(46) 



The curve falls rather abruptly to zero for {uj — u ) = 1.7uo/(nN w ) and never 
regains significant size for large frequency deviations. So, for large value of 
diffraction parameter, the frequency response curve is f(x). This is also true 
for the practically important case in which the diffraction parameter is close to 
unity. Let's examine the implication of this results for a real klystron amplifier. 
Suppose that the electron beam is in the undulator for a reasonable length, say 
for 100 periods. Then we can calculate that the FWHM of the amplification 
bandwidth is equal to (Ao;)fwhm/^o — l-l/iV w ~ 1 %. 

We have, then, a clear picture of the nature of the operation of the klystron 
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amplifier with nonzero detuning parameter, and a calculation of the band- 
width to be obtained from it. A detail analytical theory of the amplification 
at C 7^ and finite transverse electron beam size is rather complicated, but 
the discussion which we have given here is enough for ordinary purposes. 

5. 3 Output of klystron amplifier 

The amplification process in the klystron amplifier can be divided into two 
stages, linear and nonlinear. During the linear stage of amplification, signif- 
icant growth of the beam modulation amplitude and of the electromagnetic 
field amplitude takes place. Nevertheless, the beam modulation is much less 
than unity in the cascade undulators, and the largest fraction of the radia- 
tion power is produced in the output undulator, when the beam modulation 
becomes about unity. 

Suppose we have a modulated electron beam at the exit of the last cascade. 
The problem of electromagnetic wave radiation in the output undulator refers 
to a class of self-consistent problem. There is considerable analogy between 
the nonlinear mode operation of the conventional FEL amplifier and klystron 
amplifier. In the high gain linear regime of a conventional FEL amplifier the 
typical length of the change in the radiation field amplitude is about the gain 
length l g ~ r _1 , where for the one dimensional model 



_2c-f 2 z -fI A o- 2 _ 

In the case of a uniform output undulator, the bunched beam effectively in- 
teracts with the electromagnetic wave along a length which is of the order 
of the gain length l s . At this stage of the self- interaction electrons lose the 
visual fraction of their energy which results in the violation of the resonance 
condition. As a result, the beam is overmodulated, most electrons fall into 
the accelerating phase of the ponderomotive potential and the electron beam 
starts to absorb power from the electromagnetic wave. This is analogous to 
the situation with a conventional FEL amplifier. Like the efficiency of the 
conventional FEL amplifier, the efficiency r\ of the klystron amplifier with a 
long enough uniform output undulator is limited by the value of the efficiency 
parameter p: 

V-P = ci 2 J/u , 

It should be note that the saturation efficiency of the klystron amplifier does 
not depend on the amplitude of the input signal when the klystron amplifier 
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operates in the high-gain regime, i.e. when a- ni /a out <C 1. 

For some purposes it is convenient to express p in a different form 

P 



' C7 2 " 


2/3 




1/3 


" K 2 ' 


1/3 


2ua_ 








1 + K 2 _ 





(47) 



Thus the efficiency of the klystron amplifier, defined in this particular way, 
depends on only the one-third power of peak current, a very weak dependence. 
Let us make a calculation of r] for some cases. Suppose Iq = 100 A, e n = 2npm, 
7 = 1000, A w = 3 cm, K = 1.42, focusing beta function is equal to 3 m; then by 
equation (47) it follows that rj ~ 0.13 %. Returning to (47), we could imaging 
the injector linac running at a Iq — 2.5 kA. As we mentioned above, the 
dependence of 77 on the peak current is rather weak. As a result, the efficiency 
is increased by only a factor three (77 ~ 0.4 %) for Iq = 2.5 kA. 

Let us present a specific numerical example for the case of a fourth-generation 
light source. One of the key parameters to compare different radiation sources 
is their brilliance, which is simply given by the spectral flux divided by the 
transverse photon phase space. We illustrate with numerical example the po- 
tential of the proposed klystron amplifier scheme for the SASE FEL at the 
TESLA Test Facility accelerator. The average brilliance of the conventional 
TTF SASE FEL at a wavelength of A ~ 30 nm is about 

5 av ~ 10 22 photons/s/0.1%BW/mm 2 /mrad 2 



for the case when J = 2.5 kA [1]. Average photon flux and average brilliance 
varies simply as efficiency rj. Since _B av oc 77, the dependence of B &v on the 
electron peak current is contained completely in the term p. If I = 100 A, ap- 
propriate substitution in (47) shows that the average brilliance of the cascade 
klystron amplifier without bunch compression at a wavelength of A ~ 30 nm 
is about 



73 av ~ 3 • 10 21 photons/s/0.1%BW/ mm /mrad . 



On the other hand, at a wavelength of 30 nm brilliance of a third-generation 
synchrotron light source is equal to 

73 av ~ 10 17 photons/s/0.1%BW/ mm /mrad . 



A comparison to synchrotron light sources shows the drastic improvement of 
average brilliance at the klystron sources. The average brilliance of the klystron 
facility operating without bunch compression in the injector linac surpasses the 
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spontaneous undulator radiation from third-generation synchrotron radiation 
facilities by four or more orders of magnitude. However, it is not only the high 
average brilliance which makes a klystron amplifier facility a unique research 
tool in the VUV and soft X-ray regime, the radiation from this facility shows 
a very high peak brilliance too. Average brilliance is normalized to seconds at 
the highest possible repetition rate while the peak brilliance is the brilliance 
scaled to the length of a single pulse. Another type of current dependence are 
expected for peak brilliance, a 4/3 power dependence -B pC ak oc I^ 3 . The peak 
brilliance of the conventional TTF SASE FEL at a wavelength of A ~ 30 nm 
is about 

5 peaj£ ~ 10 29 photons/s/0.1%BW/ mm /mrad . 

On the other hand, the peak brilliance of the cascade klystron amplifier with- 
out bunch compression at a wavelength of A ~ 30 nm is about 

5peak - 10 27 photons/s/0.1%BW/ mm /mrad . 

Decreasing the peak current also decreases the peak brilliance of the SASE 
FEL radiation by about of factor 100, but this is still 6 orders of magnitude 
higher than that of 3rd generation synchrotron radiation sources. 

A cascade klystron amplifier can be modified to increase significantly the spec- 
tral brightness of output radiation. A reliable method to increase the klystron 
amplifier efficiency consists in an adiabatic change of the output undulator 
parameters (or, in other words, by the use of so-called undulator tapering). 
When the electron bunch passes through the last dispersion section the energy 
modulation leads to an effective (nonlinear) compression of the particles. Then 
the bunched electron beam enters a tapered output undulator, and from the 
very beginning produces strong radiation because of the large spatial bunch- 
ing. The strong radiation field produces a ponderomotive well which is deep 
enough to trap the particles, since the original beam is relatively cold. The 
radiation produced by these captured particles increases the depth of the pon- 
deromotive well, and they are effectively decelerated. As a result, much higher 
power can be achieved than for the case of a uniform undulator. 

In [14] we have described the tapering undulator scheme for a SASE FEL. 
Simulations using the code FAST provide a "full physics" description of the 
process. Despite the original spiking seeding the process of the electron density 
modulation, we effectively trap a significant fraction of the particles, and can 
achieve much higher power than for the case of an untapered output undulator. 
Another important feature of the radiation from a tapered undulator is the 
significant suppression of the sideband growth in the nonlinear regime. This 
means that in the proposed scheme the spectral brightness of the radiation 
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is increased proportionally to the radiation power. In the case of a uniform 
undulator the peak brightness is reached at the saturation point and is then 
reduced due to sideband growth. 



6 Klystron amplifier start-up from shot noise 

Till now we have considered the electron beam as a continuous medium when 
describing the theory of klystron amplifier. To some extent this is idealization, 
since in reality the beam current is produced by a large number of moving elec- 
trons. If we consider the microstructure of the electron current, we find that 
electrons enter the undulator randomly in time and space. So, we can expect 
that the FEL amplifier should possess intrinsic noise properties. Here it is rel- 
evant to remember that FELs form a separate class of vacuum-tube devices. 
The analysis of the noise properties of traditional vacuum-tube amplifiers has 
always been an important problem. This has been mainly connected with the 
practical need for reducing the intrinsic noise of the amplifier. The result of 
the experience obtained during these investigations shows that there always 
exists the fundamental effect of shot noise originating from the random emis- 
sion of the electrons from the cathode. When we analyze this effect for the 
parameters of traditional microwave amplifiers, we find that it is complicated. 
In particular, suppression of the shot noise in some frequency band can occur 
due to space charge effects. Besides the shot noise effect, there are a num- 
ber of different sources of noise which influence the operation of traditional 
vacuum-tube amplifiers. 

Fluctuations of the electron beam current serve as the input signal in the 
XFEL. As for the FEL amplifier operating in short wavelength range, its 
noise properties are defined only by the shot noise. At optical frequencies the 
quieting effect of space-charge limitation seems negligible. An FEL amplifier 
which starts up from shot noise is frequently known as a self-amplified spon- 
taneous emission (SASE) FEL. However, it is worth mentioning that such an 
essentially quantum terminology does not reflect the actual physics of the pro- 
cess. The amplification process in the SASE FEL has its origin in the density 
fluctuations in the electron beam. The latter effect is completely classical. 

The shot noise is the pure fluctuations in number which correspond to the 
fact that photoemission is a random process. The rate at which electrons are 
emitted from a photocathode is not constant in time. The emission process 
is a random one and it is impossible to predict the time dependence of the 
instantaneous current. Any such random fluctuations in the beam current 
correspond to an intensity modulation of the beam current at all frequencies 
simultaneously - including, of course, the frequency to which the undulator is 
tuned. When the electron beam enters the first undulator, the presence of the 
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beam modulation at frequencies close to the resonance frequency initiates the 
process of radiation. 

The emission of electrons from the cathode is believed to be a Poisson pro- 
cess, and from this assumption alone the total fluctuation in current can be 
deduced. Why is this the right rule, what is the fundamental reason for it, 
and how is it connected to anything else? The explanation in deep down in 
quantum mechanics. When electromagnetic fields are incident on a photosur- 
face, a complex set of events can transpire. The major steps in this process 
can be identified as absorption of a quantum light energy (i.e. a photon) and 
the transfer of that energy to an excited electron, transport of the excited 
electron to the surface, and finally, release of the electron from the surface. 
When laser light having a deterministic variation of intensity over space and 
time is incident on a photosurface, the fluctuations of the photons obey Pois- 
son statistics. This explains the relation between fundamental photon "shot 
noise" and shot noise in the electron beam. 

When describing the physical principles, it is always important to find a model 
which provides the possibility of an analytical description without loss of es- 
sential information about the features of the random process. In the first anal- 
ysis of the problem, we adopt some rather simplifying assumptions that are 
only occasionally met in practice. We will investigate the klystron amplifier 
start-up from shot noise in the framework of the one-dimensional model which 
assumes the input shot noise and output radiation have a full transverse co- 
herence. This assumption allow us to assume that the input shot noise signal 
is defined by the value of total beam current. In reality the fluctuations of 
the electron beam current density are uncorrelated in the transverse dimen- 
sion. Using the notion of the beam radiated modes, we can say that many 
transverse radiation modes are radiated when the electron beam enters the 
undulator. The one-dimensional model can be used for the calculations of sta- 
tistical properties of transversely coherent output current. In practice such an 
assumption is valid for the current at the exit of klystron with a thin electron 
beam. During the amplification process in a klystron, the number of trans- 
verse modes decreases, and the contribution of the coherent radiation to the 
total radiation power is increased up to full coherence. With this assumption, 
attention can be concentrated completely on temporal coherence effects. 

We study the case when the initial modulation of the electron beam is defined 
by the shot noise and has a white spectrum. Since in the linear regime all the 
harmonics are amplified independently, we can use the results of the steady- 
state theory for each harmonic and calculated the corresponding Fourier har- 
monics of the output current. The electric current in the time domain, I(t), 
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and its Fourier transform, are connected 

1 °° 

I(t) = — / I{uj) exp(-iotf) dw . 



The klystron amplifier has a response which can be economically expressed as 
a frequency response curve G(u). Thus, the Fourier harmonic of the current 
at the dispersion section exit and the Fourier harmonic of the input current 
are connected by the relation: 

/outM = G(u)I(u) uj > . (48) 



When uj < the Fourier harmonic is defined by the relation I*(ou) = 

It is convenient to isolate explicitly the slowly varying complex amplitude 

I(t) - Jo = I(t)e- iuot + C.C. , (49) 



where 



00 

I{t)e' iuJat = — ( I(iu)e- iujt duJ . 
2tt J 




The real output current consists of a carrier modulation of frequency ujq sub- 
jected to random amplitude and phase modulation. The Fourier decomposition 
of the output current contains frequencies spread about u . It is not possible in 
practice to resolve the oscillations in I(t) that occur at frequency of the carrier 
modulation. A good experimental resolving time is of order of picoseconds, at 
least four orders of magnitude too long to detect oscillations at the frequency 
ujq. It is therefore appropriate for comparison with experiment to average the 
theoretical results over a cycle of oscillation of the carrier modulation. The 
deviation 51 = I(t) — I has a zero cycle-average. According to (49), the cycle 
average of the square of the current deviation is 

(5I) 2 = 2 I /(*) | 2 . (50) 



Shot noise is a random statistical process and statements about shot noise are 
probability statements. These are most easily handled using the concept of 
a statistical ensemble, drawn from statistical mechanics. If the physical con- 
struction of the injector and set of the external parameters (e.g supply voltage, 
photoinjector laser beam parameters) which specify its state are known, then 
these parameters define a definite statistical ensemble to which the injector 
belongs. This ensemble consists of identical copies of the injector subject to 
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identical macroscopic external conditions. Any external property of the injec- 
tor, e.g. its laser light on the photosurface, will depend, however, not only 
on these given parameters but also on the exact microscopic internal state 
of the photocathode, and each copy of the injector in the ensemble will lead 
to a different output. We shall always denote ensemble averages by (•••). In 
practice the XFEL injector is required to produce a single electron bunch at 
a repetition rate of 100 - 100 000 bunch/s and averaging symbol (• • •) simply 
means the ensemble average over electron bunches. 



6.1 Shot noise in a linear amplifier 

The output fluctuations, expressed as current fluctuations, have a mean square 
value 

oo oo 

((5I ont ) 2 } = ^JdujJ dco'G(Lo)G*(co')(I(uj)r(uj')} . 





To calculate correlation of the spectral components of the output beam cur- 
rent, we should consider a microscopic picture of the electron beam current 
at the undulator entrance. We start the analysis for the case of a rectangular 
electron pulse of finite duration T and then we go over to the limit of an in- 
finitely long pulse. The electron beam current is made up of moving electrons 
randomly arriving at the entrance of the undulator: 

= HO E 

fc=i 

where S(- ■ •) is the delta function, (— e) is the charge of the electron, N is the 
number of electrons in a bunch and tk is the random arrival time of the electron 
at the undulator entrance. The beam current averaged over an ensemble of 
the bunches can be written in the form: 

m)= i - Z ^ L = -Io for ~\<t<\. 

The electron beam current I(t) and its Fourier transform I {to) are connected 
by: 

oo N 



46 



1 00 N 



Now we can calculate the first-order correlation of the complex Fourier har- 
monics and I(ou'): 

INN \ 

(I(cu)I(cu')) = e 2 {J2J2 exp(i^, - \u't n ) ) . 

\fc=ln=l / 

Expanding this relation, we can write: 
(/»/>')) 

= e2 ( H exp[i(^ - u')t k \ \ + e 2 ( exp(icjt fc - \uj't n ) \ 

\k=l I \k^n I 

N 

= e 2 ^2(exp[i(u - uj')t k }} + e 2 ^ (exp(i wt fc ))(exp(- icj't„)) . (51) 

k=l kj^n 

The probability of the arrival of an electron during a time interval (t, t + dt) 
is equal to dt/T. It is easy to find that 



1 J z 

(exp(i^ fe )) = ^ f e^dt k = F{u), (52) 



T/2 

I 

-T/2 

where 

=, . fuT\~ l fuoT 
F{.)={-) ,n(- 

Substituting (52) to the (51), we obtain: 

= e 2 NF(uj - u') + e 2 N(N - l)F(w)F*(w') . (53) 

When 

N | F(u) | 2 < 1 , (54) 

we can write the following expression for the first-order spectral correlation: 
(J(cu)/V)> = e 2 NF(u - u') . (55) 
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Let us consider the case of an infinitely long electron pulse with the homoge- 
neous linear density and current Iq. In this limit (T — > oo, N — > oo, eN/T = 
const.) the first-order correlation of the complex Fourier harmonics has the 
form: 

(J(cu)fV)) = e/ lim [TF(cj - = 27rel 5(uj - to') . (56) 



The following representation of the delta function has been used here: 
r/ \ sin[(wo - w)T ^ 

S(uj -u) = — 7-^ f— as T -> oo . (57) 

n(u - u) 



The definition given in (57) can be used to prove the basic property of the 
delta function. Let f{ui) be any function of u that is non-singular at uj = u , 
and consider the integral 

/ f(v)5{v -w)dw = -hm / /(w) — ? r — dw . 

J 71 T^oo J (UJq — UJ) 



Changing the variable of integration to x = (w - ou )T gives 



(u>2— u>o)T 



/ i f sm x 
f(uj)5(uj Q -uj)&uj = - lim / f((x/T)+u ) da: 

«>i (wi-wo)T 

oo 

1 f sin x 

= -f(u ) / da; = f(u ) , proved ui < u < uj 2 ■ (5£ 

7T J X 



Inclusion of the delta function in an integral therefore picks out the value of 
the integrand at the point specified by the constant in the delta function. 

Thus, if a current of mean value Iq flows in the input of an amplifier whose 
frequency response curve is G(u), the output fluctuations are 

oo 

((5I out ) 2 ) = 2(| 7(0 | 2 ) = 1° / | G{u) | 2 dco . (59) 

7T J 




We can also write (59) in the form 
el 



((Wout) 2 ) = — I G(cj ) | 2 Au A , (60) 

71 
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where the effective amplification bandwidth of the klystron is defined as 



We can further express this result by saying that the mean square value of the 
output is equivalent to a fluctuating current of mean square value 



at the input. Note that it has not been necessary to make any assumptions 
about the amplification process. The result is completely general; that is, it 
applies for any type of linear vacuum-tube amplifier [15]. 

With the preceding results in hand, it should now be possible to estimate, 
for example, the amplitude of initial beam density modulation a in . For this 
purpose it is convenient to express the result (59) in a different form. Noting 
that ((5L m ) 2 } = (af n )io and Ac<ja/^o — we can write 



where N c = N\N w /2 is the number of cooperating electrons, and N\ = 
2nIo/(euj ) is the number of electrons per radiation wavelength. If I = 100 A, 
A = 30 nm, and N w = 100, appropriate substitutions in (61) show that the 
mean square value of input is about (af n ) ~ 3 • 10~ 7 . 

Now we return to the question about the electron bunch profile. Equation 
(53) is not only true at the arbitrary chosen rectangular bunch profile, but of 
course it is true at any other bunch profile, and thus under the condition (54) 
in the long electron bunch limit the output intensity fluctuations are all the 
same. (The final statement (59) does not involve the bunch profile function, 
which appeared only in the intermediate arguments.) Let us discuss the region 
of validity of the approximation (54) for different bunch profile functions. If 
we consider the case of rectangular bunch profile, we find that the region of 
applicability of condition (54) is less than that for the case of a Gaussian 
bunch of the same duration. This is due to the fact that the bunch form factor 
| F(u) | 2 , decreas more slowly with the increase in frequency. On the other 
hand, in a realistic situation there is no sharp boundary of the bunch and the 
beam current falls to zero during some time interval A<r T <C T. When the 
beam current at the edge falls in accordance with a Gaussian law, Aot must 
obey the following conditions: 




{(8I in ) 2 ) = eI Acu A /(n) 



K) 



eu /(I 7iN w ) = 1/N ( 



(61) 



Aa T /T < 1 , Aa T u > 1 , iV/(Acx T w) 4 (wr) 2 < 1 . 
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When toT ~ 10 4 and uAa T ~ 10 2 the value (A(x t uj) 4 {ujT) 2 is equal to 10 16 . 
As a rule, the number of particles in the bunch, N, is less than or about 10 11 , 
so condition (54) is fulfilled. 

6.2 Temporal coherence 

The shot noise in the electron beam is a Gaussian random process. The 
klystron amplifier operated in linear regime can be considered as a linear filter 
which does not change the statistics of the signal. As a result, we can define 
general properties of the beam density modulation after the dispersion section 
without any calculations. For instance, the square of instantaneous amplitude 
of density modulation fluctuates in accordance with the negative exponential 
distribution. We can also state that the spectral density of the electron beam 
modulation and first-order time correlation function should form a Fourier 
transform pair (this is the Wiener Khintchine theorem). 

The correlation between the output currents at times t and t' has the form: 



(5I out (t)5I out (t')) = (/out(*)Ct(0> exp[- icu (t - 0] + CC. 



Using (48) and (56) we calculate the correlation between the complex ampli- 
tude I(t) and /*(*'): 

(J out (t)/: ut (t / ))exp[-i^ (t-t / )] 

oo oo 

= — fdoo I dec/ exp(- iut + iu't')G(u)G*(u')(I(u)I*(u')) 












el 



oc 



/ 



dccexp[— iu>(t — t')] | G(w) 



(62) 



o 



We define the first-order time correlation function as follows: 




(I(t)P(t')) 




Using (50), (59) and (62) we can write 



oo 



fdu\ G{uj) | 2 exp[-i(w-w )(t-O] 



9i(t-t') 



o 



OC' 



fdu\G(u) | 2 
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The integral should go from to oo, but is so far from ujq that the curve is all 
finished by that time , so we go instead to minus oo - it makes no difference. 
If we do that we get 

oo 

/ d(Acu) I G(Alo) I 2 exp[-i(Aw)(t - t')\ 

9At-t') = - 55 . (63) 

/ d(Acu) | G(Alu) |2 

— CO 

where Au = (cu — u ). Therefore, the slowly varying correlation function and 
the normalized spectrum of the narrow band-signal are a Fourier transform 
pair. It is seen from latter expression that the first-order correlation function 
possess the property gi(t — t') = g\(t' — t). When the klystron gain curve 
is symmetrical with respect to the resonance frequency, u , the function g 1 is 
real. In our approximation, the klystron gain function | G(Au) | 2 is symmetric 
with respect to Alu and the integration therefore needs to be performed only 
from to oo with the result being doubled. Inserting (46) into the integrand 
(63) we get 

oo 

/ f(x) cos(xr/a)dx 

9i(r) = 55 . 

Sf(x)dx 



Here we use the following notation a = irN w /uio, t — t — t'. The integral 
is a Fourier transform of f(x). In Fig. 7 we present a plot of first-order time 
correlation function. First, let us notice a remarkable feature of that plot. 
The correlation function is zero for | r | > 2a = N w X/c. Of course we could 
predict such a result. Remembering that the synchronism takes place when 
the electromagnetic wave advance the electron beam by one wavelength at one 
undulator period (see section 3), we see that 2a is the total relative slippage 
of the wave with respect to the electron beam. 

In many applications it is desirable to have a precise and definite meaning 
for the term "coherence time." Such a definition can be made in terms of 
the complex correlation function, but there are a multitude of definitions in 
terms of g\ (t — t') that can be imagined. However, in future discussions there 
is one definition that arises most naturally and most frequently. Following the 
approach of Mandel, we define the coherence time, r c , as 

oo 

t c = J | gi(r) | 2 dr . 

— oo 

Its magnitude is of the order of the inverse of the output frequency spread. 
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Fig. 7. First-order time correlation function 
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Fig. 8. Time dependence of the cycle- averaged output intensity modulation corre- 
sponds to the shot noise fluctuations in the input beam current, obtained from a 
computer simulation. The dashed line shows the mean square value of the output 
current fluctuations averaged over a time long compared to r c 

In all theory that follows, attention is restricted to modulated electron beams 
whose frequency spreads are small compared with the mean frequency, that 
is, where ujqT c is very much larger than unity. Figure 8 illustrates the type of 
fluctuations that occur in the cycle-averaged intensity modulation. The figure 
has been constructed by a computer simulation of the first cascade output in 
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which the summation is carried out explicitly for the real number of electrons 
in the electron bunch. The instantaneous intensity modulation varies in an 
unpredictable manner about the average value, and it is this kind of variation 
that characterizes the phenomenon as noise. It is seen that substantial changes 
in the square value of amplitude occur over a time span r c , but this quantity 
is reasonably constant over time intervals At <C r c . The value of the coherence 
time for the case of klystron is given by r c = 0.94a = 0A7N w X/c. With the 
numerical values N w = 100 and A = 30 nm the value of coherence time is 
r c ~ 5 fs. 

Our discussion of the current fluctuations would be incomplete if we did not 
refer the reader to another method of average calculation. The klystron output 
coherence time, r c , is in femtosecond range and I = const, over the time span 
r c is a good assumption in all practical problems. Clearly we can assume 
that the current in question is adequately modelled as a ergodic and hence 
stationary random process. A stationary process obeys statistical laws and 
is subject to conditions that do not vary with time. For a stationary process 
time and ensemble averages are equivalent. In this case averaging symbol (■ • •) 
simply means the average over a period much longer compare with r c . 



7 Further applications of klystron 

In the preceding sections we have given a simple discussion of a soft X-ray 
cascade klystron amplifier. We have not only neglected many features of the 
klystron amplifier, but we have also completely omitted some important ap- 
plications. We shall mention three here: the frequency multiplier, a scheme for 
femtosecond experiments, and multi-user soft X-ray facility. 

7. 1 Klystron frequency multipliers 

The remaining application of the klystron which we shall mention is the fre- 
quency multiplier. We have seen, in section 3, that the bunched beam at large 
values of the bunching parameter has not only a fundamental radiation fre- 
quency component, but also has considerable intensity in its harmonics. It is 
then possible to have an input undulator operating at one frequency, and an 
output undulator operating at a multiple of this frequency. The radiation in 
the output undulator will then be excited by the harmonic component in the 
electron beam, and the klystron will operate as a combination of frequency 
multiplier and amplifier. A possible application would be to the production of 
very high frequencies. A klystron gets more and more difficult to construct, as 
the frequency becomes higher and higher, on account of the reduced period of 
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the undulators. If we use a frequency multiplier klystron to reach the highest 
frequencies, it is only the output undulator which must be made with the 
smallest period; the cascade undulators, operating at a lower frequency, can 
be made with larger period, and this can accommodate an injector with larger 
emittance and make the source easier to construct. 

The idea of using two undulators, with the second undulator resonant to one 
of the harmonics of the first one (the so-called "after-burner" scheme), was 
considered in [16-18]. The technique of using an conventional FEL amplifier 
as a frequency multiplier relies on the fact that spatial bunching in the ex- 
ponential gain regime can be very strong. The first uniform undulator is long 
enough to reach saturation and produce strong spatial bunching in harmon- 
ics. The bunched beam generates coherent radiation in the second undulator 
which immediately follows the first one. The main problem with this approach 
is the large induced (correlated) energy spread which will be generated by the 
bunching of the electron beam at the fundamental frequency. While this en- 
ergy spread is necessary for the bunching, it degrades the performance of the 
radiation section at the harmonic frequency. 

The problem to be solved is how to prolong the interaction of the bunched 
electron beam with the electromagnetic wave in the harmonic radiator undu- 
lator. The present section is concerned with frequency doubler, designed for 
conditions approximating maximum efficiency A reliable method to increase 
the frequency doubler efficiency is to use a cascade klystron in which the input 
undulator operates at the fundamental resonance, the penultimate undulator 
operates at the fundamental harmonic, and the output undulator operates at 
the second harmonic. In particular it can be shown that a so-called 1-1-2 SASE 
scheme may operate with rather high efficiency in the soft X-ray wavelength 
range. (The numbers 1,2 here designate the resonant harmonic number for 
each undulator.) The klystron gain is controlled in such a way that the max- 
imum energy modulation of the electron beam at the penultimate undulator 
exit is about equal to the local energy spread. This is much smaller than the 
induced energy spread of the electrons at the exit of a saturated conventional 
SASE FEL. As a result, the operation of a frequency doubler is based on 
prolonged interaction of the bunched electron beam with the electromagnetic 
wave in the harmonic radiator. In contrast to the after-burner method, in a 
two-cascade 1-1-2 scheme frequency doubler operation is based on prolonged 
interaction of the bunched electron beam with the electromagnetic wave in 
the output undulator and copious coherent emission at double frequency is 
produced. The doubler efficiency can be optimized by monitoring the output 
radiation intensity at double frequency as a function of the compaction factor 
of the first and second dispersion section. In [14] we have described an effective 
frequency doubler for SASE FELs. For the first time the frequency multipli- 
cation scheme has been analyzed for SASE FELs. Computer modeling with 
time-dependent code FAST has demonstrated that the final output power of 
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a 1-1-2 klystron doubler is no smaller than that which could be produced by 
a 2-2-2 klystron with the same electron beam parameters. 



7.2 Scheme for femtosecond experiments 

Time-resolved experiments are used to monitor time-dependent phenomena. 
The standard technique for high-resolution time-resolved experiments is the 
pump-probe scheme in which a process is started by a short pulse of radia- 
tion (pump) and the evolution of the process is then observed (probed) after 
a delay by means of a second pulse of radiation, generally at another pho- 
ton energy. The obvious temporal limitation of the visible pump/X-ray probe 
technique is the duration of the X-ray probe. Here we will concentrate on the 
performance which can be obtained by a klystron XFEL operating without 
bunch compression in the injector linac. At these sources, the X-ray pulse du- 
ration is about 10 ps. This is longer than the timescale of many interesting 
physical phenomena. The new principle of pump-probe techniques described 
below offers a way around this difficulty 

Our studies have shown that the soft X-ray cascade klystron holds great 
promise as a source of radiation for generating high power single femtosec- 
ond pulses. Our femtosecond soft X-ray facility concept is based on the use of 
an X-ray SASE FEL combined with a femtosecond quantum laser [19]. The 
operation of a femtosecond soft X-ray SASE FEL is illustrated in Fig. 9. An 
ultrashort laser pulse is used to modulate the density of electrons within the 
femtosecond slice of the electron bunch at a frequency uj opt . The seed laser 
pulse will be timed to overlap with the central area of the electron bunch. 
This ultrashort laser pulse serves as a seed for a modulator which consists 
of an uniform undulator and a dispersion section. The interaction of seed 
pulse with the electron beam produces an energy modulation at cj opt . This en- 
ergy modulation is converted into spatial bunching in the dispersion section. 
Density modulation at the modulator exit is about 20 %. The energy modu- 
lation introduced by the modulator is smaller than the local energy spread. 
Following the modulator the beam enters the first undulator of the klystron 
amplifier which is resonant with soft X-ray radiation at frequency c^o- The 
process of amplification in the soft X-ray klystron amplifier develops in the 
same way as was explained in previous sections: fluctuations of the electron 
beam current density serve as the input signal. At the chosen level of density 
(energy) modulation the SASE process develops nearly in the same way as 
with a non-modulated electron beam because of the large ratio of cooperation 
length to optical wavelength. By the time the beam is bunched in the klystron 
cascades, at frequency uq, the density modulation amplitude has reached sat- 
uration. This leads to amplitude modulation of the density at the sidebands 
ujq ± u opt . The sideband density modulation takes place at the part of the 
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Fig. 9. Scheme for time-resolved experiments based on the generation of femtosecond 
pulses by sideband-seeded soft X-ray cascade klystron amplifier. In this scheme, a 
femtosecond soft X-ray pulse is naturally synchronized with the femtosecond optical 
pulse from the seed laser and cancels jitter 

electron pulse denned by the duration of the seed laser pulse, which is much 
shorter than the electron pulse. Following the last cascade the beam enters 
the output undulator which is resonant at the frequency u — oo opt . Because 
the beam has a large component of bunching at the sideband, coherent emis- 
sion is copiously produced within the femtosecond slice of the electron bunch. 
Analytical methods are of limited use in to the study of the sideband-seeded 
SASE FEL and numerical simulations must be used. Simulations with time- 
dependent code FAST provide a "full physics" description of the process in a 
sideband-seeded SASE FEL. The results of numerical simulations confirm our 
simple physical considerations [19]. 

We assume that the SASE bandwidth is much less than the separation of 
the sidebands from the main peak. This requirement is of critical importance 
to the overall performance of the femtosecond SASE FEL facility. In this 
case, the (output) undulator can be used to distinguish the sideband density 
modulation from the density modulation at the central frequency. Obviously, 
this requirement is easier to achieve for long radiation wavelengths. For 300 nm 
laser radiation and for 30 nm output radiation, for example, the amplification 
bandwidth of the cascade klystron must be much smaller than 10 %. This 
condition may be easily satisfied in practice. The long (unbunched) electron 
pulse is one of the advantages of the adopted FEL design. Since the electron 
pulse duration (about 10 ps) is much longer than the time jitter of the electron 
and seed laser pulses (about ±1 ps), the synchronization of the optical laser 
with the electron pulses is not a problem. 

Pump-probe techniques which are commonly used with optical lasers, are 
highly desirable in order to make full use of the femtosecond soft X-ray pulses. 
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density modulation of the electron bunch 
at the exit of second cascade has the form 

j3cos(cjoptt + c))cos(w t + d) 




GJ 



Fig. 10. Description of the sideband generation for the case of the density modulation 
as initial condition 

Since, in this case, precise timing is needed with a jitter of less than 10 fs, in 
[19] we suggested to combine the femtosecond soft X-ray pulses generated by 
the sideband-seeded SASE FEL with optical pulses generated by the seed laser 
system. It should be emphasized that in the proposed scheme a femtosecond 
soft X-ray pulse is naturally synchronized with his femtosecond optical pulse 
and cancels jitter. This development of XFEL based pump-probe experiments 
allow us to investigate phenomena at timescale down to 100 fs. 

7.3 Multi-user distribution system for XFEL laboratory 

An X-ray laboratory should serve several, may be up to ten experimental 
stations which can be operated independently according to the needs of the 
user community. On the other hand, the multi-user distribution system has 
to satisfy an additional requirement. Passing the electron bunch through the 
bending magnets must avoid emittance dilution due to coherent synchrotron 
radiation (CSR) effects. For very short bunches and very high peak current, 
CSR can generated energy spread in the bending magnets and thus dilute 
the horizontal emittance. As a result, the prefered layout of a conventional 
SASE FEL is a linear arrangement in which the injector , accelerator, bunch 
compressors and undulators are nearly collinear, and in which the electron 
beam does not change direction between accelerator and undulators. 

The situation is quite different for the klystron amplifier scheme proposed 
in our paper. Since it operates without bunch compression in the injector 
linac, the problem of emittance dilution in the bending magnets does not ex- 
ist. An electron beam distribution system based on unbunched electron beam 
can provide efficient ways to generate a multi-user facility - very similar to 
present day synchrotron radiation facilities. A possible layout of a soft X-ray 
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Fig. 11. Diagram of a possible fourth- generation synchrotron facility using sponta- 
neous emission free-electron laser klystron amplifiers 

FEL laboratory based on an electron ring distribution system is shown in Fig. 
11. The layout of the laboratory follows a similar approach as that used for 
synchrotron light sources. The SASE FEL user facility consists of an injector 
linac and electron beam distribution system. The injector is composed of a 
RF gun with photocathode and a main (superconducting) linac. In order to 
make efficient use of the new source it is proposed to segment the full cir- 
cumference of a distribution system into arcs which are repeated a number of 
times to form a complete ring. Each cell includes a two-cascade klystron and 
a bending magnet. A specific realization of the electron ring cell is sketched in 
Fig. 12. The electron beam transport line guiding electrons from the injector 
linac to the experimental hall is connected tangentially to one of the straight 
sections of the electron ring. In order to obtain a useful separation between 
the experimental areas behind the photon beam lines, an angle of 36 degrees 
between two neighboring lines would be desirable. Thus, ten beam lines can 
be installed on a complete electron ring. Using klystron (electromagnetic) dis- 
persion sections in each cell as switching elements it is possible to quickly 
switch the FEL photon beam from one experiment to the other, thus provid- 
ing multi-user capability. Users can define the radiation wavelength for their 
experiment independently of each other to a very large extent, since they use 
different undulators. Injector linac and electron beam transport lines operate 
at fixed parameters. At a fixed electron energy the magnet gap of the klystron 
undulators can be varied mechanically for wavelength tuning. This design 
makes it possible to make various wavelengths of SASE radiation available in 
the XFEL laboratory quasi- simultaneously. It is a great advantage that injec- 
tor and electron beam transport lines in the new scheme of multi-user facility 
operate at fixed parameters and that an "electron switchard" is not required. 

When a relativistic electron beam passes through the undulator, it emits in- 
coherent radiation. This process leads to an increase of the energy spread in 
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Fig. 12. Electron ring cell design. Using a klystron electromagnetic dispersion sec- 
tion as a switching element it is possible to quickly switch off (on) the cell klystron 
amplifier thus providing multi-user capability. This design makes it possible to 
make various wavelengths of SASE radiation available in the XFEL laboratory 
quasi-simultaneously 

the beam due to quantum fluctuations of the undulator radiation. This effect 
grows significantly with an increase in electron beam energy, strength field 
and length of the undulator. This effect should be carefully taken into ac- 
count when designing multi-user distribution system. For the case of planar 
undulator the expression for the rate of energy diffusion can be written in the 
following form [20]: 



dz 15 2tt w V ; ' 



(64) 



where A c is the Compton wavelength, r e is the classical electron radius, F(K) ~ 
1.2K for K ^> 1. The next numerical example illustrates the amplitude of the 
effect. If 7 = 10 3 , K = 1.42, A w = 3 cm, and the total undulator length is 
equal to 100 m, appropriate substitution in (64) shows that the energy spread 

increase due to quantum fluctuations is about \J ((AS) 2 ) ~ 1 keV at the end 
of the 10th klystron amplifier, and has a negligible effect on the klystron per- 
formance. 



8 Conclusion 



The high-gain klystron amplifier described in this paper is an attractive alter- 
native to other FEL configurations for operation in short wavelength range. 
An distinguishing feature of the klystron amplifier is the absence of apparent 
limitations which would prevent operation without bunch compression in the 
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injector linac. We have illustrated the proposed cascade klystron scheme us- 
ing the parameters of the TESLA Test Facility. Although the present work 
is concerned primarily for use in the soft X-ray spectrum, its applicability is 
not restricted to that range, for example an X-ray SASE FEL is a suitable 
candidate for application of cascade klystron scheme described here. 
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